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Abstract 

A dynamical coupled-channel model is presented for investigating the nucleon resonances (N*) 
in the meson production reactions induced by pions and photons. Our objective is to extract the 
N* parameters and to investigate the meson production reaction mechanisms for mapping out the 
quark-gluon substructure of N* from the data. The model is based on an energy-independent 
Hamiltonian which is derived from a set of Lagrangians by using a unitary transformation method. 
The constructed model Hamiltonian consists of (a) Ty for describing the vertex interactions 
N* <-> MB,7tttN with MB = ^N,i:N,r]N,TTA, pN,aN, and p <-> tttt and a <-» tttt, (b) V22 for 



(N 
> 



the non-resonant MB — * M'B' and nir — * tttt interactions, (3) vmb,ttttN for the non-resonant 
MB — > 7T7T./V transitions, and (4) v^N^N for the non-resonant ttttN — ► ttttN interactions. By 
applying the projection operator techniques, we derive a set of coupled-channel equations which 
satisfy the unitarity conditions within the channel space spanned by the considered two-particle 
MB states and the three-particle ttttN state. The resulting amplitudes are written as a sum of 
non-resonant and resonant amplitudes such that the meson cloud effects on the ./V* decay can be 
explicitly calculated for interpreting the extracted N* parameters in terms of hadron structure 
calculations. We present and explain in detail a numerical method based on a spline- function ex- 
pansion for solving the resulting coupled-channel equations which contain logarithmically divergent 
one-particle-exchange driving terms Z MB M , B , resulted from the ttttN unitarity cut. This method 
is convenient, and perhaps more practical and accurate than the commonly employed methods of 
contour rotation/deformation, for calculating the two-pion production observables. For complete- 
ness in explaining our numerical procedures, we also present explicitly the formula for efficient 
calculations of a very large number of partial-wave matrix elements which are the input to the 

coupled-channel equations. Results for two pion photo-production are presented to illustrate the 

( p) 

dynamical consequence of the one-particle-exchange driving term Z MB M , B , of the coupled-channel 
equations. We show that this mechanism, which contains the effects due to ttttN unitarity cut, can 
generate rapidly varying structure in the reaction amplitudes associated with the unstable particle 
channels 7rA, pN, and aN, in agreement with the analysis of Aaron and Amado [Phys. Rev. D13, 
2581 (1976)]. It also has large effects in determining the two-pion production cross sections. Our 
results indicate that cautions must be taken to interpret the N* parameters extracted from using 
models which do not include 7T7rA cut effects. Strategies for performing a complete dynamical 
coupled-channel analysis of all of available data of meson photo-production and electro-production 
are discussed. 

PACS numbers: 13.60.Le, 13.60.-r, 14.20.Gk 
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I. INTRODUCTION 



With the very intense experimental efforts at Jefferson Laboratory (JLab), Mainz, Bonn, 
GRAAL, and Spring-8, extensive data of photo-production and electro-production of ir, 
r], K, u, (p, and two pions have now become available 0]. Many approaches have been 
developed accordingly to investigate how the excitations of nucleon resonances (N*) can be 
identified from these data. The objective is to extract the A"* parameters for investigating the 
dynamical structure of Quantum Chromodynamics (QCD) in the non-perturbative region. 
The outstanding questions which can be addressed are, for example, how the spontaneously 
broken chiral symmetry is realized, and how the constituent quarks emerge as effective 
degrees of freedom and how they are confined. In this work, we are similarly motivated and 
have developed a dynamical coupled-channel model for analyzing these data. 

The irN and 7 A" reaction data in the A^* region are most often analyzed by using two 
different kinds of approaches. The first kind is to apply the models which are mainly the 
continuations and/or extensions of the earlier works. These include the analyses by us- 
ing the Virginia Polytechnic Institute-George Washington University (VPI-GWU) Model 
(SAID)0, the Carnegie- Mellon-Berkeley (CMB) modelQ, and the Kent State University 
(KSU) model 4]. Apart from imposing the unitarity condition, these models are very phe- 
nomenological in treating the reaction mechanisms. In particular, they assume that the 
non-resonant amplitudes, which are often comparable to or even much larger than the reso- 
nant amplitudes, can be parameterized in terms of separable or polynomial forms in fitting 
the data. Furthermore, their isobar model parameterizations do not fully account for the an- 
alytical properties due to the nnN unitarity condition, as discussed, for example, by Aaron 
and Amado[5|. We will address this important question later in this paper. 

The second kind of analyses account for the reasonably understood meson-exchange mech- 
anisms. For numerical simplicity in solving the scattering equations, they however neglect 
the off-shell multiple-scattering dynamics which determines the meson-baryon scattering 
wavefunctions in the short range region where we want to map out the quark-gluon substruc- 
ture of N*. The inrN unitarity condition is also not satisfied rigorously in these analyses. 
The most well-developed along this line are the Unitary Isobar Models (UIM) developed 
by the Mainz group (MAID)p and the Jlab-Yeveran collaboration 0], K-matrix coupled- 
channel models developed by the Giessen groupjH and KVI group0, and the JLab-Moscow 
State University (MSU) model of two-pion production. More details of these approaches 
have been reviewed recently in Ref.0|. 

As we have learned recently in the A region, the results from the approaches described 
above are useful, but certainly not sufficient for making real progress in understanding the 
structure of N* states. For example, the empirical values of A^-A transitions extracted by 
using SAID and MAID are understood within the constituent quark model only when the 
very large pion cloud effects are identified in the analyses based on dynamical models(10, 
30. The essence of a dynamical model is to separate the reaction mechanisms from 



the internal structure of hadrons in interpreting the data. To make similar progress in 
investigating the higher mass A^*, it is highly desirable to extend such a dynamical approach 
to analyze the meson production data up to the energy with invariant mass W ~ 2 GeV. This 
is the objective of this work. Our goal is not only to extract the resonance parameters, but 
also to interpret them in terms of the current hadron structure calculations. The achievable 
goal at the present time is to test the predictions from various QCD-based models of baryon 
structure. It is also important to make connections with Lattice QCD calculations. The 
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FIG. 1: The total cross section data of meson production in jp reaction. Left: 1 — ir and 2 — it 
production are compared. Right: KY ( K + A, K + TF > , K°T, + ), r/p, and ujp production are compared 
with some of the 1 — tt and 2 — ir production 



Lattice QCD calculations are now being carried out [13] to give a deeper understanding of the 
iV-A transition. A systematic Lattice QCD program on N* is also under development [13] . 

The main challenge of developing dynamical reaction models of meson production reac- 
tions in the N* region can be seen in FigGJ We see that two-pion photo-production cross 
sections shown in the left-hand-side become larger than the one-pion photo-production as 
the 7p invariant mass exceeds W ~ 1.4 GeV. In the right-hand-side, KY ( K + A, K + T,°, 
K°T, + ), rjp, and up production cross sections are a factor of about 10 weaker than the 
dominant 7T + 7T~p production. From the unitarity condition, we have for any single meson 
production process 'jN -> MB with MB = ttN, r]N, uN, KA, KY, 

i( T iN,MB - T* MBnN ) = T lN ,M' B' Phi' B'T* MB M , B , 

M'B' 

J t-T^N,TTTTNPlTTTNT MBt7T7TN , (l) 

where p a denotes an appropriate phase space factor for the channel a. The large two-pion 
production cross sections seen in FigC] indicate that the second term in the right-hand- 
side of Eq.(£Q) is significant and hence the single meson production reactions above the A 
region must be influenced strongly by the coupling with the two-pion channels. Similarly, 
the two-pion production 7^ — > 7r7riV is also influenced by the transition to two-body MB 
channel 

i(T^N,n-KN - T* nN ^ N ) = T lN,M'B'PM'B'T* nNyMIBI 

M'B' 

(2) 

Clearly, a sound dynamical reaction model must be able to describe the two pion production 
and to account for the above unitarity conditions. 

The development of meson-baryon reaction models including two-pion production channel 
has a long history. It was already recognized in 1960 's, as discussed by Blankenbecler and 
Sugar [l5l| . that the dispersion-relation approach, which has been v ery successful in analyzing 
the data of ttN elastic scattering[l6| and 7iV — > ttN reactions ["LA Il8| . can not be used 
to analyze the data of two-pion production. The reason is that apart from the tiN and 
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irirN unitarity cuts, it is rather impossible to even guess the analytic structure of two- 
pion production amplitudes. Furthermore, the dispersion relation models are difficult to 
solve because of their bi-linear structure which is the price of only dealing with the on-shell 
amplitudes. 

Ideally, one would like to find alternatives to analyze the ttN and reaction data 
completely within the framework of relativistic quantum field theory. The Bethe-Salpeter 
(BS) equation has been taken historically as the starting point of such an ambitious approach. 
The complications involved in solving the BS equation have been known for long time. For 
example, its singularity structure and the associated numerical problems were very well 



discussed in Refs.[19, 20j, |21|. The BS equation contains serious singularities arising from 
the pinching of the integration over the time component. In addition to the two-body 
unitarity cut, it has a selected set of n-body unitarity cuts, as explained in great detail in 



Ref . [21[ . Considerable numerical efforts are already needed to solve the Ladder BS equation 
for nN elastic scattering, as can be seen in the work of Lahiff and Afnanji^]. Using the Wick 
rotation, they can solve the Ladder BS equation below two-pion production threshold with 
very restricted choices of form factors. It is not clear how to extend their work to higher 
energies. 

The first main progress in finding an alternative to the dispersion-relation approach 
was perhaps also made by Blankenbecler and Sugar [l5|. By imposing the unitarity con- 
dition, they show that the Bethe-Salpeter equation can be reduced into a covariant three- 
dimensional equation which is linear and can be managed in practice. Compared with the 
dispersion relation approach, the challenge here is account for the off-shell dynamics. This 
approach was later further developed by Aaron, Amado, and Young (AAY)[23(. With the 
assumption that all interactions are due to the formation and decay of isobars, they de- 
veloped a set of covariant three-dimensional equations for describing both the 7rN elastic 
scattering and ttN — > nnN reaction. They however had only obtained 23[ 24], HL Ig| a very 
qualitative description of the nN data and only investigated very briefly the electromagnetic 
meson production reactions. Their results suggested the limitation of the isobar model and 
the need of additional mechanisms. For example, the N* excitation mechanisms are not 
included in their formulation. They then proposed an approach to include the additional 
mechanisms phenomenologically in fitting the data by using the "minimal" equations which 
are rigorously constrained by the ttN and itkN unitarity conditions and have the correct an- 
alyticity of the isobar model. The AAY approach was later applied mainly in the studies of 
iiNN systems, such as those by Afnan and Thomas [271] and by Matsuyama and Yazaki Jiil] . 
Development in this direction was well reviewed in Ref. |29| . 

The dynamical study of ttN scattering was pursued further in 1980's by Pearce and 



Afnan [30L l3lL |32j. They derived the tcN scattering equations by using a diagrammatic 
method, originally developed for investigating the nNN problem[29], to sum the perturba- 
tion diagrams which are selected by imposing the unitarity condition. Furthermore, they 
relate the 7rN scattering to the cloudy bag model by extending the work of Thieberg, Thomas 
and Miller (33l. 34l 35l | to include the ttttN unitarity condition. 



Since 1990 the ttN and 'yN reactions have been investigated mainly by using ei- 
ther the three-dimensional reductions [36] of the Bethe-Salpeter equation or the unitary 
transformation methods [Tol. l3?| . These efforts were motivated mainly by the success of 
the meson-exchange models of NN scattering 3^] , and have yielded the meson-exchange 



models developed by Pearce and Jennings |39j], National Taiwan University- Argonne Na- 
tional Laboratory (NTU-ANL) collaboration ji^, 41 1, Gross and Surya[42], Sato and 
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Lee(To|, [U, Julich GroupEM E IH Fuda and nis collaborators 53, and Utretch- 
Ohio collaboration |48l l49j |. All of these dynamical models can describe well the data in 
the A region, but have not been fully developed in the higher mass N* region. The main 
challenge is to include correctly the coupling with the ttttN channels. 

We now return to discussing the two-pion production channel which is an essential part 
of our formulation. Most of the recent two-pion production calculations are the extensions 
of the isobar model of Luke and S6ding[50]. The production mechanisms are calculated from 
tree-diagrams of appropriately chosen Lagrangians. The calculations of Valencia Group |5l| 
included the tree diagrams calculated from Lagrangians with 7, N, 71", p, A(1232), iV*(1440), 
and iV*(1520) fields. To describe the total cross section data in all charged mrN channels, 
they also included |H3| the production of A(1700) and p effect aris ing from A/"* (1520). 

The model developed by Ochi, Hirata, Katagiri, and Takaki p3L IHil 155^ contains the 
tree diagrams calculated from Lagrangians with 7, tt, p, u>, N, A and iV*(1520) fields. 
An important feature of this model is to describe the excitation of iV*(1520) within an 
isobar model with three channels ttN, pN, and it A. They found that the invariant mass 
distributions of all charged channels of 'jp — > mrN can be better described if the pseudo- 
scalar ttNN coupling is used. They also found that the iV*(1520) — > pN decay is the 
essential mechanism to explain the differences between the invariant mass distributions of 
7T + 7r° and 7r°7r°. Similar tree-diagram calculations of two pion photo-production have also 
been performed by Murphy and LagetjH^]. 

The analyses |5 71 l58l . |59( of two pion production by using the JLab- Moscow State Uni- 
versity (JLAB-MSU) isobar model considered only the minimum set of the tree diagrams 
proposed in the original work of Luke and S6ding[5(j. However, they made two improve- 
ments. They included all 3-star and 4-star resonances listed by the Particle Data Group 
and used the absorptive model developed by Gottfried and Jackson 6(| to account for the 
initial and final state interactions. They found that the ttNA form factor is needed to get 
agreement with the data of jp — > 7r~A ++ , while the initial and final state interactions are 
not so large. In analyzing the two-pion electro-production data, they further included a 
7nrN phase-space term with its magnitude adjusted to fit the data. This term was later 
replaced by a phenomenological particle-exchange amplitude which improves significantly 
the fits to the data. With this model, they had identified 5^| a new iV*(| + , 1720) and the 
production of the isobar channel 7r + .Di3(1520) which has never been considered before. 

The common feature of all of the two-pion production calculations described above is that 
the coupled-channel effects due to the unitarity condition, such as that given in Eqs.(JT])-(|2J), 
are not included. The problems arise from this simplification were very well studied by 
Aaron and Amado^, and will be discussed later in this paper. While the results from these 
tree-diagram models are very useful for identifying the reaction mechanisms, their findings 
concerning N* properties must be further examined. To make progress, it is necessary to 
develop a coupled-channel formulation within which the nuN channel is explicitly included. 
In this paper, we report our effort in this direction. 

We have developed a dynamical coupled-channel model by extending the model developed 



in Refs. |l0l lllj to include the higher mass N* and all relevant reaction channels seen in Fig^ 



Our presentations will only include two-particle channels MB = jN, ttN, rjN and three- 
particle channel miN which has resonant components 7rA, pN, and aN. But the formulation 
can be easily extended to include other two-particle channels such as uN, KA and K~E and 
three-particle channels such as m]N and KKN . 

Our main purpose here is to give a complete and detailed presentation of our model and 
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the numerical methods needed to solve the resulting coupled-channel equations. A complete 
coupled-channel analysis requires a simultaneous fit to all of the meson production data 
from iiN and 7 A reactions, such as the total cross section data illustrated in FigJT] and 
the very extensive data from recent high precision experiments on photo-production and 
electro-production reactions. Obviously, this is a rather complex problem which can not 
be accomplished in this paper. Instead, we will apply our approach only to address the 
theoretical questions concerning the effects due to 7r7riV unitarity cuts. For this very limited 
purpose, we present results from our first calculations of 7 A — ► ttttN reactions. 

In section II, we present the model Hamiltonian of our formulation. It is derived from a 
set of Lagrangians, given explicitly in Appendix A, by applying the unitary transformation 
method which was explained in detail in Refs. (loL l6l| . The coup led- channel equations are 
then derived from the model Hamiltonian in section III with details explained in Appendix 
B. In section IV, we explain the procedures for performing numerical calculations within 
our formulation. The numerical methods for solving the coupled-channel equations with 
7nrN cut are explained in section V. Results of 7p — ► ttttN are presented and discussed in 
section VI. A summary and the plans for future developments are given in section VII. For 
the completeness in explaining our numerical procedures, several appendices are given to 
present explicitly the formula for efficient calculations of a very large number of partial- wave 
matrix elements which are the input to the coupled-channel equations, and to explain how 
the constructed resonant amplitudes are related to the information listed by the Particle 
Data Group (PDG)|fi3. 



II. MODEL HAMILTONIAN 

In this section we present a model Hamiltonian for constructing a coupled-channel re- 
action model with 7A, ttN, rjN and 7T7rA channels. Since significant parts of the tihN 
production are known experimentally to be through the unstable states 7rA, pN, and per- 
haps also aN, we will also include bare A, p and a degrees of freedom in our formulation. 
Furthermore, we introduce bare N* states to represent the quark-core components of the 
nucleon resonances. The model is expected to be valid up to W = 2 GeV below which three 
pion production is very weak. 

Similar to the model of Refs. 1~D| (commonly called the SL model), our starting point is 



a set of Lagrangians describing the interactions between mesons (M = 7, tt, rj, p,u,a- ■■) and 
baryons (B = N,A,N* ■ ••). These Lagrangian are constrained by various well-established 
symmetry properties, such as the invariance under isospin, parity, and gauge transformation. 
The chiral symmetry is also implemented as much as we can. The considered Lagrangians 
are given in Appendix A. By applying the standard canonical quantization, we obtain a 
Hamiltonian of the following form 

H = J h(x, t = 0)dx 

= H + Hj, (3) 

where h(x, t) is the Hamiltonian density constructed from the starting Lagrangians and the 
conjugate momentum field operators. In Eq.([5|). i^o is the free Hamiltonian and 

Hi — Fmb^b' + X! h,M'M"+^M , (4) 

M,B,B' M.M'.M" 
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where Tmb^b' describes the absorption and emission of a meson(M) by a baryon(B) such 
as ttN <-> N and ttN <-> A, and h,M'M"*-*M describes the vertex interactions between mesons 
such as 7T7T <-> p and 77r <-> 7r. Clearly, it is a non-trivial many body problem to calculate 
meson-baryon scattering and meson production reaction amplitudes from the Hamiltonian 
defined by Eqs.(j2J)-(l3J)- To obtain a manageable reaction model, we apply a unitary trans- 
formation method [loLl6lj| to derive an effective Hamiltonian from Eqs.([5|)-P|). The essential 
idea of the employed unitary transformation method is to eliminate the unphysical vertex 
interactions MB B' with masses mjf + ms < itib 1 from the Hamiltonian and absorb 
their effects into MB — > M'B' two-body interactions. The resulting effective Hamiltonian is 
energy independent and hence is easy to be used in developing reaction models and perform- 
ing many-particle calculations. The details of this method have been explained in section II 
and the appendix of Ref.|lO|. 

Our main step is to derive from Eqs.p])-(|1}) an effective Hamiltonian which contains inter- 
actions involving tttcN three-particle states. This is accomplished by applying the unitary 
transformation method up to the third order in interaction Hj of Eq.p|). The resulting 
effective Hamiltonian is of the following form 

H eff = H + V, (5) 

with 

= (6) 



where K a = \jm? a + p* a 2 is the free energy operator of particle a with a mass m a: and the 
interaction Hamiltonian is 

V = T v + v 22 + v , (7) 

where 

IV = {XXX! Tn*~,MB + IV— ttttAt) + ^M*-*7r7r} + {c.C.} , (8) 

N* MB M* 

V22 = X V MB,M'B' + Vtttt ■ (9) 
MB, M'B' 

Here c.c. denotes the complex conjugate of the terms on its left-hand-side. In the above 
equations, MB = jN, irN, rjN, it A, pN, aN represent the considered meson-baryon states. 
The resonance associated with the bare baryon state N* is induced by the vertex interactions 
^n*-*mb and r jy* — >7rnN • Similarly, the bare meson states M* = p, a can develop into 
resonances through the vertex interaction hM*-+wir- These vertex interactions are illustrated 
in Fig|2fa). Note that the masses Mjy* and mP M , of the bare states iV* and M* are the 
parameters of the model which will be determined by fitting the irN and tttt scattering data. 
They differ from the empirically determined resonance positions by mass shifts which are 
due to the coupling of the bare states with the meson-baryon scattering states. It is thus 
reasonable to speculate that these bare masses can be identified with the mass spectrum 
predicted by the hadron structure calculations which do not account for the meson-baryon 
continuum scattering states, such as the calculations based on the constituent quark models 
which do not have meson-exchange quark-quark interactions. It is however much more 
difficult, but more interesting, to relate these bare masses to the current Lattice QCD 
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(b) v 22 = jxx^ + 

B ^» B' it t 



(c) v 23 = 




FIG. 2: Basic mechanisms of the Model Hamiltonian defined in Eqs.(8)-(10). 
M M' 

**">-:"'' :M n 



B B n B' 



v s v u v c 



FIG. 3: Mechanisms for vmb,m'B' of Eq.(9): (a) direct s-channel, (b) crossed u-channel, (c) one- 
particle-exchange t-channel, (d) contact interactions. 

calculations which can not account for the scattering states rigorously mainly because of the 
limitation of the lattice spacing. 

In Eq.flOJ), Vmb,m'B' is the non-resonant meson-baryon interaction and v niT is the non- 
resonant 7T7T interaction. They are illustrated in FigJ2fb). The third term in Eq.l jTj) describes 
the non-resonant interactions involving ttitN states 

v' = V 2 3 + V 33 (10) 

with 

V 23 = ^[{vmB.ihtn) + (CC.)] 
MB 

v 33 — V-k-kN^-kN ■ 

They are illustrated in FigJ2fc). All of these interactions are defined by the tree-diagrams 
generated from the considered Lagrangians. They are illustrated in Fig|3] for two-body in- 
teractions vmb,m'B' and in FigEJfor Vmb^-kN- Some leading mechanisms ofv^ and v^n^n 
are illustrated in FigJHl The calculations of the matrix elements of these interactions will 
be discussed later in the section on our calculations and detailed in appendices. Here we 
only mention that the matrix elements of these interactions are calculated from the usual 
Feynman amplitudes with their time components in the propagators of intermediate states 
defined by the three momenta of the initial and final states, as specified by the unitary 
transformation methods. Thus they are independent of the collision energy E. 
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FIG. 4: Examples of non-resonant mechanisms of vmn,-k-kN with M = ir or 7 (denoted by lon£ 
dashed lines). Mj denotes the intermediate mesons (ir,p, u). 
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FIG. 5: Examples of non-resonant mechanisms of v nn and v^N-KirN 



III. DYNAMICAL COUPLED-CHANNEL EQUATIONS 



With the Hamiltonian defined by Eqs.(jo])-()10p , we follow the formulation of Ref. 63|| to 
define the scattering S-matrix as 

S ab (E) = 5 ab -(2ir)iT ab (E), (11) 

where the scattering T-matrix is defined by 

T ab (E) =< a\T(E)\b> 

with 

T(E) = V + V l,. T(E) . (12) 
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FIG. 6: Graphical representations of Eqs.(14)-(19). 



Since the interaction V, defined by Eqs.fl7j)- (fTU|) . is energy independent, it is rather straight- 
forward to follow the formal scattering theory given in Ref.|63j] to show that Eq. (J12|) leads 
to the following unitarity condition 

{T{E)-T\E)) ab = -2m^Tl(E)5(E c -E)T cb (E), (13) 

c 

where a, b, c are the reaction channels in the considered energy region. 

Our task is to derive from Eq. (fT2j) a set of dynamical coupled-channel equations for 
practical calculations within the model space iV* © MB © ttttN. In the derivations, the 
unitarity condition Eq.(|13|) must be maintained exactly. We achieve this rather complex 
task by applying the standard projection operator techniques jH^], similar to that employed 
in a study|65( of ttNN scattering. The details of our derivations are given in Appendix B. To 
explain our coupled-channel equations, it is sufficient to present the formula obtained from 
setting TN*->mrN = in our derivations. The resulting model is defined by Eqs.(B74)-(B96) 
of appendix B. Here we explain these equations and discuss their dynamical content. 

The resulting MB — > M'B' amplitude Tmb-*m'b> in each partial wave is illustrated in 
FigEl It can be written as 

TmB,M'B'{E) = tMB,M'B> M'B' (E) , (14) 

The second term in the right-hand-side of Eq. ()14|) is the resonant term defined by 

t MBM'B'( E )= F MB ^ N *(E)[D(E)]i !j r N j-+M>B'(E) , (15) 

N*,N* 
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with 



[D(E)-%(E) = (E — M° N: )6 h3 - £ M -(£) , (16) 
where M^* is the bare mass of the resonant state N*, and the self-energies are 

= X! Fn*->mbGmb{E)Tmb^n*{E) . (17) 

MB 1 3 

The dressed vertex interactions in Eq.(fT5j) and Eq.(fI7|l are (defining Tmb^n* = ^n*^mb) 

f MB-+N* (E) = VmB-^N* + X! ^MB,M'B'{E)Gm'B i {E)Ym'B'-^N* , (18) 

M'B' 

Fn*->Mb(E) = Tn*^MB + X! ^N*->M'B'GAI'B'(E)tM'B',MB(E) . (19) 

M'B' 

The meson-baryon propagator Gmb in the above equations takes the following form 

Gmb{e) = ~f k — k 1 y TjaTT, ' (20) 

t, - K B - K M - ^MB{ti) + xe 

where the mass shift Y*mb{E) depends on the considered MB channel. It is Hmb{E) = for 
the stable particle channels MB = ttN, r]N. For channels containing an unstable particle, 
such as MB = it A, pN, aN, we have 

Y>mb{E) = [< MB\g v - y , . g\?\MB >] un - connected (21) 

with 

9V = ^A-^TrN + hp^lT + ^<7-*7T7T • (22) 

In Eq. (l21|) "tm — connected" means that the stable particle, tt or N, of the MB state is a 
spectator in the ttttN propagation. Thus Timb(E) is just the mass renormalization of the 
unstable particle in the MB state. 

It is important to note that the resonant amplitude t M , B , MB {E) is influenced by the 
non-resonant amplitude tu'B' ,mb(E), as seen in Eq. (fTB]) - (fT^j) . In particular, Eqs.(fT% j) - (fT§|) 
describe the meson cloud effects on N* decays, as illustrated in Fig|7|for the A — > decay 
interpreted in Refs. [ToL 11 1. This feature of our formulation is essential in interpreting the 



extracted resonance parameters. 

Here we note that the N* propagator D(E) defined by Eq. ffTtjj) can be diagonalized to 
write the resonant term Eq. (fTo^) as 

+R (U JJ\ _ V" ^MB^N*( k ) T N^M'B'(k') . s 

t MB, M'B' Wi-L; E - Mft.(E) + iTf.(E) ' {2S) 

where Tmb^n* an d mass parameters M^* (E) and T^i (E) are of course related the dressed 
vertexes T N *_> MB and self energies Sjj defined in Eqs.flr7 |) - (fT^|) . Eq.(J22I) is similar to the 
usual Breit-Wigner form and hence can be used to relate our model to the empirical resonant 
parameters listed by Particle Data Group. This non-trivial subject is being investigated in 
Ref. 1661. 
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pion cloud effect 



FIG. 7: Graphical representation of the dressed Fa^n interpreted in Refs. H], . 

The non-resonant amplitudes t MB>M i B i in Eq. (fT^|) and Eqs. (fTHj) - (fT^j) are defined by the 
following coupled-channel equations 

tMB,M',B' (E) = V MB ,M'B' 

(E)+ V MB ,M"B" \E)Gm"B" \E)tM"B",M'B' (E) (24) 

M"B" 

with 

Vmb,m'b' (E) — vmb,m>b> + Zmb,m>b> (E) . (25) 

Here Zmb,m'B'(E) contains the effects due to the coupling with 7T7riV states. It has the 
following form 

Z M b,m'b>(E) = <MB\ F- — ^ —Ft | M'B' > 

E — H — v nnN + te 

[3 MB, M'B' ^-*MB 

(E)) (26) 

with 

(27) 

F = gv + vmb,ttttN , (28) 

where gv has been defined in Eq. (l22j) . Note that the second term in Eq. (J26|) is the effect 
which is already included in the mass shifts £mb of the propagator Eq.(j20I) an d must be 
removed to avoid double counting. 

The appearance of the projection operator P n7T N in Eqs. (}2~Tj) and (}2l)j) is the consequence of 
the unitarity condition Eq. ffTHj) . To isolate the effects entirely due to the vertex interaction 
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z (E) = ^\ 

/cn fF^ TP^ 

FIG. 8: One-particle-exchange interactions Z^ nA (E), Z pNnA and Z^ NwA of Eq.(30). 

9v — Ta^ttAt + hp^-K + — »7r7r ; we use the operator relation Eq.(B33) of Appendix B to 
decompose the ttttN propagator of Eq. (|26J) to write 

Zmb,m'b'(E) = Zm B)M , b ,{E) + Zm B}M , b ,(E) . (29) 

The first term is 

Zmb,M'b>{E) = < MB | g v P ™ N , g\ \ M'B' > -[5 MB)M >b^ mb (E)} . (30) 

Obviously, Z MB M , B , (E) is the one-particle-exchange interaction between unstable particle 
channels 7rA, pN, and aN, as illustrated in FigJBl The second term of Eq. ()29j) is 

Z$ B , M 'B'(E) = < MB | F — W^(E) Ft | M'B' > 

+ < MB | | B'M' > 

+ <MB\ VMB N —Jjgi—gi v | M'B' > 
E — H + ie 

+ < MB | ^mb.^tv „ . v{ IB ^ N | M'B' > . (31) 

E — H + ie 

Some of the leading terms of Z$ B M , B , (E) are illustrated in FigEl Here t n7r N,imN(E) is a 
three-body scattering amplitude defined by 

tinvN,TVKN{E) = VttttN + V^JSf— — — — ; — —V-jnrN (32) 

E - K n - K v - K N - v^ N + ie 

where v n7T jy has been defined in Eq. (127)1 . Few leading terms of Eq.([32|) due to the direct 
s-channel interaction v s (illustrated in Fig. 3) of v-kn^n are shown in FigllUI These terms 
involve the nnN propagator 1/(E — — K n — Kn + ie) and obviously can generate mr N 
cut effects which are due to the ttttN vertex. This observation indicates that the irN scat- 



tering equation of Aaron, Amado, and Young 23j can be related to our formulation if the 
interactions which are only determined by the nnN vertex are kept in the equations pre- 
sented above. We however will not discuss this issue in this paper. The relations between our 
formulation and the AAY model can be better understood in our next publication j6?J where 
we will determine the strong interaction parts of our Hamiltonian by fitting ttN reaction 
data up to invariant mass W = 2 GeV. 

The amplitudes Tmb,m'B' = ^mb,m'B' + ^mb m'b> defined by Eq. (fT3)l can be used directly 
to calculate the cross sections of ttN — > nN, r)N and jN — > ttN, rjN reactions. They are 
also the input to the calculations of the two-pion production amplitudes. The two-pion 
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N 



+ 



+ 




+ 

FIG. 9: Examples of mechanisms included in Z$ B M i B r(E) of Eq.(31). 




+ 

FIG. 10: Some of the leading order terms of t W7r N,irnN of Eq.(32). The open circle represents the 
direct s-channel interaction v s illustrated in Fig. 3 for the MB = M'B' = ttN case. 

production amplitudes resulted from our derivations given in Appendix B are illustrated in 
FiglTTl They can be cast exactly into the following form 

T nwN ,MB(E) = T%" NMB (E) + T^ NMB (E) + T^ NMB (E) + T^ NMB (E) (33) 
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7iN,7utN 'MB.tcN jcN.tcjcN 




T, 



MB,7iA 



T, 



MB,p(a)N 



FIG. 11: Graphical representations of T n7r N t MB defined by Eqs.(33)-(37). 



pdir 

- ttttN,MB 



(E) 



T. 



with 

< ^iJv^)! ^2 V ttttN,M'B'[^M'B',MB 
M'B' 

+GM'B'(E)(tM'B' ,mb(E) + t^ IIB , MB )]\MB > , 

< 4-J N (E)\ri^ N G^ A (E)[^ A , MB (E) +t* A>MB (E)]\MB > 

T£ NMB {E) = <4~MmU^GMm^,MB(E)+t^ NMB (E)]\MB> 
In the above equations, the tttiN scattering wave function is defined by 

where the scattering operator is defined by 



sfc&(E) =< mrN\[l+t 



tvttN.tttvN 



(E) 



l 



E - K„ - K„ - K N + ie 



(34) 
(35) 

(36) 

(37) 

(38) 
(39) 



Here the three-body scattering amplitude t n7rN ^ nN (E) is determined by the non-resonant 



interactions u™., v. 



j u -kN,-kN 



and v „ 



tN.ttttN 



as defined by Eq.flUZ). 



We note here that the direct production amplitude T%£ N MB (E) of Eq. fjMj) is due to 
VnnN,MB interaction illustrated in FigEl while the other three terms are through the un- 
stable 7rA, pN, and aN states. Each term has the contributions from the non-resonant 
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amplitude tu'B 1 ,mb(E) and resonant term t^, B , MB (E). As seen in Eq.(fT5 |) -(fT!? j) . the reso- 
nant amplitude t^ rB , MB [E) is influenced by the non-resonant amplitude tM>B',MB(E). This 
an important consequence of unitarity condition Eq.()13|). 



IV. CALCULATIONS 

The N* information can be accurately extracted only when the extensive meson produc- 
tion data of ttN and 7A reactions are analyzed simultaneously. Obviously, this is a rather 
complex task by using the dynamical coupled- channel formulation described in section III. 
In addition, it is a highly non-trivial numerical task to solve the coupled-channel equa- 
tion Eq. plf which contains a logarithmically divergent driving term Zmb,m'B'(E) defined 
by Eqs. (j29j) - (f3"T|) . As a first step, we focus in this work on the development of numerical 
methods for solving this coupled-channel equation. This then allows us to perform two-pion 
photo-production calculations to investigate the effects due to the imN cut effects which are 
not included in the recent two-pion production calculations, as briefly reviewed in section I. 

To proceed, we first note that the matrix elements of Z^ BM , B ,, as defined by Eq. fpIT]) . 

is expected to be weaker than the other driving terms vmb,m'B' and Z MB M , B , because it 
involves more intermediate states. For our present purpose of developing numerical methods, 
this rather complex term can be neglected in solving the coupled-channel Eq. (j24J) . For 
simplicity, we also neglect the non-resonant interactions on the final ttttN state by setting 
< ip^~J N (E)\ ^< 7T7rA| in the calculation of two-pion production amplitudes defined by 
Eqs.(E3J-©. 

To make contact with recent experimental developments, we focus on the 7A — > ttttN 
process. Our task is therefore to develop numerical methods for solving the following equa- 
tions 

T^NaN^E) = T^ N)lN {E) + f^ N ^ N (E) + T^ NaN (E) + T°™ N ^ N (E) (40) 



with 



T-xKN^Ni.E) — < 7171 N I v inrN,yN + V W7T n,tvnG w n (E)[U Nj yN + t nN ^ N ] 1 7 A > , (41) 

T:Xn(E) = < nnNlTi^G^mUME) + t« AaN )\jN > , (42) 
f^ NaN {E) = < ^N\h\^G pN {E)[i pNaN {E) + tfx^N > , (43) 
T^ N , lN {E) = < ™N\hU^G aN {E)[t aNaN {E) + t^wW-YN > . (44) 

Here the non-resonant scattering amplitudes Imb.m'B' is obtained from solving Eq. (j24j) with 
one of its driving term Zj^ B M , B , set to zero. To the first order in electromagnetic coupling, 
the matrix elements of these non-resonant amplitudes are calculated from the following 
coupled-channel equations 

i M B,M'B'(k,k',E) = V MBtM 'B'(k,k',E) 

M"B" J 

(45) 

tMB,yN(k, q, E) = VmB,~/N (k, q) 

+ 53 / d k'i MBi M'B'(k,k',E)G MIBI (k',E)v M 'B', 1 N(k,q) (46) 



M'B' 
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with 

V M B,M'B'{k, k', E) = v M B,M'B'(k, k') + Z { MB,M'B'(k, k', E) (47) 

where MB = ttN, rjN, 7rA, pN, aN. Despite the neglect of some of the terms of the formu- 
lation presented in section III, the calculations based on the above equations are already 
far more complex than all of existing calculations of two-pion production based on the tree- 
diagram models or K-matrix coupled- channel models. This is however a necessary step 
to correctly account for the meson-baryon scattering wavefunctions in the short range re- 
gion where we want to extract and interpret the N* parameters using the data of meson 
production reactions, as discussed in section I. 

In the following subsections, we describe our numerical procedures for solving Eqs. (j45|) - 
(}4Tj) to get the non- resonant amplitudes t mb,m' b' > calculating the resonance amplitudes 
^mb m'b'i an d evaluating the two-pion production amplitudes Eqs. pUjl - PH) . 



A. Non-resonant amplitudes 

We solve Eq. (}4l)j) in the partial- wave representation. To proceed, we follow the convention 
of Goldberger and Watson jH^] to normalize the plane-wave state \k > by setting < k\k' >= 
5(k — k'). In the center of mass frame, Eq. (jll|) then leads to the following formula of the 
cross section of M(k) + B(-k) — ► M(k') +B(-k') for stable particle channels MB, M'B' = 
7JV, ttN, r)N 

<h - * 177)2 PM>B>(k')p MB (k) . 1 — £ £ | < M'B'\T(E)\MB > | 2 



m 3M' m iB m ',» r : m ' 



3m m 

(48) 



with 



< M'B'\T{E)\MB >= 

< JM^'jitf i ^M^im ' ^B^jg > ^'b^'tb |-^M'B',Mb(^'; k, E^j \j Al177>j M j iMf^iM ) 3B^jg ■> TB^rg > , 

(49) 

where [(j'm? m j M ), (*m> m i M )] an d [C/B m j s )) ( r B m r B )] are the spin-isospin quantum numbers 
of mesons and baryons, respectively. The incoming and outgoing momenta k and k' are 
defined by the collision energy E 

E = E M (k) + E B (k) = E M ,{k') + E B '{k') , (50) 

and the phase-space factor is 

kE M (k)E B (k) 

PMB(k) = 7T . (51) 

The partial-wave expansion of the scattering amplitude is defined as 

T M 'B>,MB(k', k,E)= £ £ I ^S^(^) > 7^, B , LSMB (k f , k, E) < Y%?£0)\ 

JM,TM T LS,L'S' 

(52) 
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(53) 



where the total angular vector in the spin-isospin space is defined by 

\ Y lUmJb)s(^) > = J2 \jMm jM ,i M ,m iM ;j B m jB ,T B m TB ><TM T \iMr B m iM m TB > 

all m 

x < JM\LSm L m s >< Sm s \j M j B m jM m jB > Y LrrlL (k) . 

Clearly, Eas. (J52> (j53 jl lead to 

TjJs' 'm> 'B> ,lsmb(^ 5 ^5 E) 

= J dk' J dk< Y£§*jJ s , {k')\T M >B<,MB{ti, k; E) \ Y^™ T s (k) > . 

(54) 

By also expanding the driving term VMB,M'B'(k,k',E) of Eq.pS^l into the partial- wave 
form similar to Eq. (|52jl . we then obtain a set of coupled one-dimensional integral equations 

t J us'M'B',LSMB(k ',k,E) = V^, M , B , LSMB (k',k,E) 

+ ^2 ^2 / k u2 dk n V[7 S i M i B i L n S ii M ,i B „{k , k , E) 

M"B" L"S" 

xG M »B»(k", E)t J J ls ,, M „ B ,, tLSMB (k" ', k, E) , (55) 

where the driving term is 

^L'S'M'B',LSMB(k 5 k) = v J jJ s , M i B i LSMB {k , k) + Z^g, M , B , LSMB (k ,k,E). (56) 

The above partial- wave matrix elements of the non- resonant interaction Vm'B',mb an d one- 
particle-exchange interaction Z M , B , MB (E) are given in Appendices C and E, respectively 
There the numerical methods for evaluating them are also discussed in some details; in 
particular on the use of the transformation from the helicity representation to the partial- 
wave representation. 

The propagators in Eq. (J55j) are given in Appendix B. Taking the matrix elements of 
Eqs.(B84)-(B90), we have 

E - E M (k) - E B (k) + te 
for stable particle channels MB = 7rN, r)N, and 

G M B(k,E) = — , r , r — y- . (58) 

E - E M (k) - E B [k) - S MB (fc, E) 

for unstable particle channels MB = 7rA, pN, aN with 

SW*. E) = / fa- m J^f (q))2 + pp + fc , (59) 
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( F) 

FIG. 12: Logarithmically divergent moon-shape regions of the matrix elements of Z^ A Va(^'> ^> ^0 
(solid curves) and Z pN nA (k' ,k, E) (dashed curves). 



where the vertex function /A j7r jv(g) is from Ref.[l(|, f P ,mr(<l) and f a7T7T (q) are from the isobar 
fits[68| to the im phase shifts. They are given in Eqs.(D7)-(D9) of Appendix D. 

To solve the coupled-channel integral equation Eq. ()55)1 . we note that the matrix elements 
of their particle-exchange driving terms Z^ A ^ A (k, k', E) and Z pN ^ nA (k, k , E) (FiglBJ) contain 
singularities due to the nirN cuts. This can be seen in Eq. (E5) of Appendix E which is the 
essential component of their partial- wave matrix elements Eq.(E2). Qualitatively, they are 
of the following form 



A J S„ A (L'S',LS, l,k,k')Pi(x) 



j sva^ k' E) ~ / dx- 

LSnA,LS,AK ^J_ x E -E^k)-E,{k')-E N {k + k')+ie 

r+i A J p l^{L'S',LS,l,k,k')Pi{x) 



J J — 1 



E - E n (k) - E N (k') - E n (k + k') + ie 



(62) 
(63) 



where A JT is a non-singular function, Pi{x) is the Legendre polynomial, and x = k ■ k! . 
One can easily see that these two driving terms diverge logarithmically in some momentum 
regions. For E = 1.88 GeV, they are within the moon-shape regions of Fig^J Their 
boundary curves are defined by E - E v (k) - E^k') - E N (k ± k') = for Z$ vA and by 
E — E n (k) — E]si{k') — E w (k± k') = for Z p ^ wA . In FigH2l we show the rapid change of the 

matrix element Z\ A \ A {k, k ; E) at E = 1.88 GeV and k = 300 MeV/c when the momentum 
k is varied to cross the moon-shape region. In particular, the imaginary part (dashed line) is 
non-zero only in a narrow region. The matrix elements of Zpff wA (k, k'\ E) have the similar 
singular structure. 

With the singular structure illustrated in FigHSl Eq.([55|) can not be solved by the stan- 
dard subtraction method. To get ttN — > ttN, rjN and 'jN — > 7rN, r)N on-shell scattering 
amplitudes, it is sufficient to apply the well-developed method of contour rotation to solve 
Eq. (j53j) on the complex momentum axis defined by kg = ke~ l6 with 9 > 0. However, 
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FIG. 13: Matrix elements of the one-particle-exchange term Z^ nA (k, k , E) for L = L = 1, J 
5/2, T = 1/2 at k' = 300 MeV/c and E = 1.88 GeV. 



the resulting half-off- shell transition amplitudes tMB,-yN(ke,q', E) with MB = 7rA,pN,aN, 
defined on the complex momentum kg, can not be used directly to evaluate the matrix ele- 
ments Eqs.flJT])-(JlD) for calculating the two-pion production amplitudes. Considerable effort 
is needed to find an appropriate contour integration for getting the desired matrix elements 
on the real momentum axis. The situation is similar to the calculations of deuteron breakup 
in nd or pd reactions, as well discussed in the literatures We overcome this difficulty by 
applying the spline-function method developed in the study of nNN reactions 0,E3- This 
method is explained in details in the next section. 

The solutions of Eq. (jB3j) are then used to calculate the non- resonant photo-production 
amplitudes Eq. (J46|) . Here we use the helicity-LSJ mixed-representation that the initial 'jN 
state is specified by their helicities, A 7 , Ajv, but the final MB is defined by the (LS)J angular 
momentum variables 



VMB,yN(k,q) 



Y Y Y I Y L<jMj B )s(k) > v WMB,\^X N m TN ^-> Q'E) 
JM,TM T LS A 7 Ajv 

X Y ~4~ ~ - P M,(A 7 -A iv )(0g^g»-0g) < A 7 ,Aivm Tjv |, 



(64) 



where D J m m ,{4>,d,-4>) 
Eq. pUj) then leads to 



e t(m+m')<f,^ m ,(9) with d 3 m m ,(9) being the Wigner rotation function. 



JT 



V 



JT 

LSMB,\-y\ N m 



(k,q,E) + £ £ k' 2 dk't 

LSMB,L'S'M'B' 

{k,k',E) 

M'B' US' 



xG M 'B'(k', E)v J L 7 s , M , B , x x (k', q, E) 



(65) 

The matrix elements vlg MB ^ Nmr (k, q, E) considered in our calculations are given in Ap- 
pendix F. This unconventional representation, which is convenient for calculations, can be 
related to the usual multipole expansion, as also given in appendix G. 
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B. Resonant amplitudes 



Our next step is to calculate the resonant term defined by Eq.(|15|). Here we need to 
perform calculations using bare N* — > MB vertex functions generated from some hadron 
models. Obviously, this is a non-trivial task and beyond the scope of this work. In particular, 
one needs to analyze the consistency between the employed hadron model and our reaction 
model. Instead, we use the diagonalized form Eq. (J23|) and simply make some plausible 
assumptions to calculate the resonant amplitude t MB M , B , by using the information listed by 
Particle Data Group (PDG)[62|. In the center of mass frame we write Eq.(|23J) for 7iV — > MB 
transition in the helicity-LSJ mixed-representation as 

where M^* is the resonance position. The calculations of the decay functions Tff* LSM'B'(k') 
and r^ i A 7 A Jv m TJV (<?) are explained in appendix I. They are 



r J N T *,LSMB(k) - {2 ^, r2 -j== ^J^^J^L [GUmb] fS(k, k R )(^-) L (67) 

f^x^m^te) = (2^)3/2 y^^y^^V / ^« A AL iV ]^A T (9, q^h^-x^ , (68) 

where kn and qR are defined by M^* = -Eb(&r) + Em^r) = qR + En^r)- The form factors 
are normalized such that fls(kR,kR) = 1 and g{ T (qR,qR) = 1. For simplicity, we choose 
fH(k, k R ) = (A 2 /((A; - k R f + A 2 )) 2 and g( T (q, q R ) = (A 2 /((g - q R f + A 2 )) 2 with A = 650 
MeV/c. As explained in Appendix I, the forms Eqs.()67J ) -()68J ) are chosen such that the 
coupling strength G J L T SMB is related to the partial decay width T MB (Nj T ) of the considered 
N* -> MB 

Tmb(N* jt ) = J:\G J l I mb \ 2 , (69) 

LS 

and the 7 AT — > N* helicity amplitude is related to the partial decay width by 



q R m N 8 JT 2 JT 



,2 



Eq. (|7()jl is defined in the N* rest frame and the photon momentum q is in the quantization 
z-direction. 

The total width Y^*{E) in Eq. (j66|) is parameterized by using the variables of N* — > ttN 

decay as 

" ~ A 2 




2Z/7T 



1 Lvr+4 



_(^-M 2 + A 2 J ' 

where T^t is the value given by the Particle Data Group, L n is the orbital angular momentum 
of the considered ttN state and 

_ kE N (k)E n (k) 

P{k) ~ ^ E N (k) + E^kY (?2) 

In the above equations, A;^ is the pion momentum at energy E while ko n is evaluated at 
E = Mn*. We set the form factor parameter A = 650 MeV/c. Our main results on the 
effects due to the 7tttN cut are not changed much if we vary the cutoff A in Eqs. (}6Tj) - (fTTj) . 
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C. jN — > ttttN cross sections 

Our last step is to calculate the cross sections of 7(g) + N(p) — > 7r(ki) + 7r(k 2 ) + N(p'). 
With the S-matrix defined by Eq. tfTTj) and the normalization < k\k' >= 5{k — k'), we have 

da = ^)^5 {4 \p + q-k 1 -k 2 - p')dkidk 2 dp' 



Vrel 
1 

4 



x- E E I < h,m l ,k2,m i2 ,p , m' jN m' T jT vnNyjN (E)\qX 1 ,p\ N m TN > | 2 , (73) 



where and the isospin quantum number of the outgoing two pions, m'- and 

m' Tjv are the spin-isospin quantum numbers of the outgoing nucleon. The initial 7iV state is 
specified by their helicities A 7 ,Atv and the nucleon isospin t n . With some straightforward 
derivations, the differential cross section with respect to the tttt invariant mass M n7T can be 
written in the center of mass ( p = —q and k — {k\ + k 2 ) = —p')as 

<tn r dQ k [ dn kl2 — J a — (74) 



dM n7T J J dQkd^ik 12 dM n7T 
with 

d° = {2 ^\ EN ^ V \ \ EN ^ E ^ k ^ E ^ k2 \ \k k] 

dQ k dQ kl2 dM n7T 1 M E l[ E l[ ' 12j 

1 

4 



x 7 X X I < ^1, m n? ^2, mi 2 , P' m 'j N m ' TN \TTrnN,yN(E)\q\ y , p\ N m TN > |L> 



1 



(75) 

—* —* —* — * 

where k\ and k 2 are related to the relative momentum k\ 2 and center of mass momentum k 
of the 7T7T subsystem by a Lorentz boost 

* = & + + ij)^ 1 ' (76) 

* = + i 1 ^' 12 ' " (77) 



with 



= 2E n (k 12 ) , (78) 

^(A:) = E v (ki) + E„(k 2 ) 

= y/Afa + P , (79) 

£ = E N (k) + E^{k) . (80) 



The above equations lead to 



M 2 

k 12 = \j^-ml. (82) 

(83) 
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The matrix element < kx,mi 1 ,k 2 ,mi 2 ,p'mj N m , TN \T^ nNaN (E)\qX J ,pX^m m > can be cal- 
culated from the partial-wave matrix elements of tMB,yN(E), t R MB,-yN(E), and the vertex 
interactions A — > ttN and p,a — ► ttti. As an example, the matrix element of the term 
T^ N N {E) defined by Eq.(jI3J) can be calculated from 

m jp> m i P Is JMj,TM T LS 



x <^jp,*p^;^m ijv ,rjvm Tw |y La ^ )s T (A;) > 



fp,TTw(kl2) 



E-E N (k)-E p (k)-E pN (k,E) 

-LS , pA r ,A 7 A]vm rjv 



x KW.A^ ?: *0 + ^$v,A 7 A*-nv 9, #)] (84) 



where ly/^^g'Cp) > h as been defined in Eq.fl53J), j n = = and hence only s = and 
/ = j p are allowed in the sum. 

Expressions similar to Eqs. (|75|) - (j%4"j) can be easily obtained for the differential cross sec- 
tions with respect to the irN invariant mass M^n by changing the labels of variables. 



V. NUMERICAL METHODS 

To illustrate the numerical method we have developed for solving the coupled-channel 
equation Eq lJHHJ) with a singular particle-exchange driving term Z^ B M , B ,, it is sufficient 
to consider the Alt-Grassberger-Sandhas (AGS) integral equation[72j within a simple three 
identical bosons model of Amado [73l| . This model describes the scattering of a boson b from 
a two-boson bound state d via a d — ► bb form factor g(q) = go/(q 2 + /3 2 ) with q denoting 
the relative momentum between the two outgoing bosons. The form factor is normalized as 
/ k 2 dkg 2 (k)/(B + ^-) 2 = 1 with B being the binding energy of the two-boson subsystem. 

After partial wave projection, the AGS equation in each partial-wave is 

X(j/,po,E) = Z(j/,po,E) + J p 2 dpZ(p',p,E)r(p,E)X(p,p ,E) , (85) 

where X(p,p , E) is the half-off-shell bd — > bd scattering amplitude. The one-particle ex- 
change driving term Z(p',p,E) and the propagator r{p,E) are calculated by using the 
familiar non-relativistic kinematics. In the center of mass system, they are 



Z(j/,p,E) = l -f^dxP L {x 



i f 1 go go 

(\p' + Ip\ 2 + P 2 )(\p+^\ 2 + P 2 ) 

i 



x ; ^ , (86) 

2m 2m 2m 

r-\p,E) = (E 2 (p,E) + B) 



x 



g 2 (k) 

(B + %)*(E 2 (p,E)-% + ie) 



1 - (E 2 (p, E) + B)l k 2 dk _ ■ J" L ■ - 1 , (87) 



where L is the orbital angular momentum, E 2 (p, E) = E— 3p 2 /4m, and Pl{x) is the Legendre 
polynomial with x = p' ■ p. 
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FIG. 14: Logarithmically divergent moon-shape region of the matrix elements Z(p' ,p, E) of Eq.(85) 
of the Amado Model, p, p' and E are in unit of U = 2m = 1 with E = 1. 

Besides the 2-body bound state pole at E 2 {p, E) + B = E - 3p 2 /(4m) + B = 0, the 
interaction Z(p',p, E) in the kernel of Eq. (j85|) has logarithmic singularity for energies above 
the three-particle breakup threshold. With the parameters H = 2m — 1, B = 1.5, (3 = 5, 
and the total energy E = 1, one can see from the energy denominator of Eq. (J86|) that the 
interaction Z(p',p,E) is singular in the moon-shape region of Figim Since the singularity 
depends on both p and p', it is difficult to solve the integral equation Eq. by using 
the standard subtraction methods. Although there are well-known methods of contour- 
deformation to avoid the singularity, we will solve the equation without contour-deformation 
by employing the interpolating function. Because mathematical problems of the singular 
integral equation (|85|) are well discussed in Ref . [zl| for example, we will concentrate on the 
practical numerical procedures. 

Let us choose appropriate grid points {p^ and write the unknown function X(p,p ,E) 
in terms of an interpolation function Si(p) 

X(p,p ,E) = Y,S t (p)X( Pl ,p ,E) . (88) 

i 

By inserting Eq.(|88|l) into eq.([85|). one obtains the matrix equation 

X{ Pj , Po ,E) = Z{p j ,po,E) + Y^K ji X{p i , Po ,E) , (89) 

i 

where 

K 3l = / p 2 dpZ( Pj ,p, E)r(p, E)Si(p) = / P 2 dpZ(p ]}P} E)r(p } E)S l (p) . (90) 

The integration in Eq. (j9U|) can be carried out as precisely as necessary since the interpolation 
functions Si(p) are known and the logarithmic singularity can be integrated as / dx\n(x) = 
xln(x) — x. The integration over the 2-body bound state pole of r(p, E) can be worked out 
by using the standard technique of pole subtraction. 

The choice of interpolation functions S%(p) depends on the property of the function to be 
interpolated. For example, the Lagrange interpolation polynomials are employed in Ref. |zij 
with some care near the breakup threshold. In the case of polynomial interpolation, however, 
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FIG. 15: Half-off-shell amplitude X(p,po, E) of Eq.(85) of the Amado Model, p and po are in unit 
of h = 2m = 1 with E = 1. The dot-dashed curves are from deriveing term Z(p,po, E) of Eq.(85). 

some changes in a small region may give rise to global effects. Therefore it is better to use 
the spline interpolation which depends locally on the grids points, i.e., the function Si(p) 
dominates around the grid point p^. Moreover, the spline interpolation is known to be less 
oscillating compared to the polynomial interpolation. 

Spline functions are defined in terms of piecewise polynomials which are connected 
smoothly over the whole region. Since cubic splines are mostly employed, we will explain 
it in some detail. There are several kinds of spline functions depending on the condition 
of continuity. Among them, natural splines and Hermitean splines are very useful. Their 
characteristic properties are: 

(1) natural splines: first and second derivatives are continuous at the grid points. It is 
a global spline in the sense that the function Si{p) depends on the whole grid points. It is 
known that the natural spline interpolation has a minimum curvature property. 

(2) Hermitean splines: Only first derivatives are continuous at the grid points. It is a 
local spline in the sense that the function Si{p) (pi < p < Pi+i) depends on 4 grid points 

{Pi-l, Pi, Pi+1, Pi+2}- 

Since the practical ways of calculating the spline functions Si(p) are well described in 
Ref. [zi| for natural splines and in Ref.[76[ for Hermitean splines, we will not repeat them 
here. 

The choice of spline functions certainly depends on the behavior of the solution X(p,p ). 
As is well-known, there appears a square-root singularity at the breakup threshold [zH]. More 
precisely, the amplitude X(p,po, E) goes like (ps — p) i+1 / 2 ( i is an angular momentum 
of the 2-body bound state) below the breakup threshold ps- Therefore, in the case of 
t = 0, the derivative is not continuous at pb and there appears a sharp change of the 
amplitude. The straightforward application of the spline interpolation is not suitable since 
it requires the smooth continuation. One of the ways to take into account this singular 
threshold behavior is to divide the whole region [0, oo] into two regions [0,ps] [pb, oo], and 
employ Hermitean spline interpolation in each region. It is also recommended that the 
grid points are suitably modified to account for the singularity near the breakup threshold, 

i.e., p' = \jp 2 B ~ P 2 (p < Pb) and p' = \Jp 2 — p 2 B (p > Pb)- In order to check the spline 
interpolation for the square-root singularity, it is a good exercise to fit the simple model 
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and to examine the accuracy of the interpolation. This exercise also will give some idea 
about the distribution of the grid points. 

Now we will explain how the spline function method works in a calculation of Eq.([85|) 
for the Amado model with the parameters h = 2m = 1, B = 1.5, (3 = 5, and the total 
energy E = 1. As discussed above, the interaction Z(p',p,E) given in Eq.tjHHJ) is sin- 
gular in the moon-shape region of Figll4l To choose the grid points for solving Eq. (j85|) 
with the input Eqs. (j86|) - (JHTj) . we first identify some typical momenta: the on-shell mo- 



mentum po = J4m(E + B)/3 of bd elastic scattering, the tips of the moon-shape re- 



gion on the coordinate axes p en d = \fmE, the breakup threshold ps = jAmE/3, and 




Pb = JmE/3 at which the moon-shape boundary has its maximum value from each co- 



ordinate axis. We then choose p a = 0, pt, = JmE/3 = 0.408, p c = p en d = y/mE = 



0.707, pd = Pb — yj4mE/3 = 0.816, p f = p = \jAm(E + B)/3 = 1.291, p g = p max = 20 
and p e = (pd + P/)/2 = 1.053. These momenta are chosen to make 6 regions as 

Ra = [Pa, Pb], Rb = [Pb, Pc],Rc = [Pc, Pd],Rd = [Pd, Pe\, Re = [Pe, Pf], Rf = \Pf,Pg\- 

In addition to the grid points of those typical momenta, we prepare {2,2,4,3,3,9} grid 
points in each region respectively, and thus 30 mesh points are used in solving the matrix 
equation Eq.(jHHJ) . They are distributed in equal space for R a , Rb and R e , while modi- 
fied grid points p\ are equally spaced near the breakup threshold for R c and Rd- In the 
region Rf, grid points are distributed as geometrical series with the ratio r = 1.5 ;i.e., 
p = 1.291, 1.456, 1.704, 2.075, 2.632, 3.468, 4.722, 6.602, 9.423, 13.65, 20 . 

In order to evaluate the integral Eq.ffDTH) accurately, we have employed 4-point Gauss- 
Legendre integration formula for each interval [p n , p n +i] which has no singularity. For the 
interval including the logarithmic singularity, we have changed the integration variable by 
explicitly taking account the location of the singularity as 



where p s (p n < p s < p n +i) is the singular point. The variable is changed as p = p s + t 3 and 
^i = {Ps —Pn) 1 ^, ^2 — (Pn+i —Ps) 1 ^- This manipulation explicitly removes the logarithmic 
divergence from the integrand. 

Thus, we have prepared two kinds of mesh points, i.e., one is the grid points {p^} at 
which the solution X(pi,p , E) is to be found by solving the matrix equation Eq. (JS§|) . and 
the other is to carry out the integration of Eq. (|9TI|) as precisely as required. 

The calculated amplitude X(p, po, E) for zero total angular momentum are the solid curve 
(real part) and dashed curve (imaginary part) shown in Fig ll5l which can be compared 
with the similar calculation of Ref. [t3] - The amplitude X(p,p ,E) is dimensionless and 
normalized as X(po,po) = (r]e 2tS — l)/(2i) at the on-shell point. One can see clearly the 
square-root singularity at the breakup threshold. We have also carried out the calculation 
with natural splines. Although natural splines are not suitable for the square-root singularity, 
it is practically possible to imitate the singularity by distributing many grid points around 
the breakup threshold. For example, the elastic amplitudes calculated by two different 
splines agree within the accuracy of 1%, since the on-shell point is away from the breakup 








(92) 
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threshold. In practice, both amplitudes coincides fairly well except for the small region 
around the breakup threshold. In Fig lTol we also show the contribution (dot-dashed curve) 
from the driving term Z(p,po,E) defined by Eq.(|86p. Its differences with the solid and 
dashed curves clearly show that the multiple scattering effects are very important. 

The method described above can be readily extended to solve the coupled-channel equa- 
tion Eq.(|55p. To be more specific, let us consider the case of E = 1.88GeV. As discussed 

( pi 

in the previous section, the partial- wave matrix elements of the driving terms Z^ nA and 

Z^ nA of Eq.([55j) diverge logarithmically in the moon-shape regions shown in FieTHZl To 
choose the grid points for applying the spline function expansion method, we first select the 
following momenta 
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The momentum p§ is the on-shell momentum of the ttN state. Pc{pi) corresponds 
to the momentum at which the invariant mass of the ttN (tttt) subsystem of the 
wN state is myi = m^ + m n (2m 7T ). This momentum can be considered as the 
"breakup" threshold of the unstable particle channels 7rA (pN and aN). Specifically, 
we take {p , pi, P2, P3, Ps, Pe, P7, Ps, Pmax} for 7r-spectator channel (ttN, it A), and 
{Po, Pi, P2, Pa, Ps, P&, Pi, Ps, Pmax} for iV-spectator channel (pN,aN). For example, 
numerical values at E — 1.88 GeV are : p\ = 80.29, p 2 = 334.8, p 3 = 448.7, p 4 = 
539.1, p 5 = 605.8, p 6 = 619.4, p 7 = 669.7, p 8 = 696.3 and p max = 6000. For 8 regions 
Ri = [Po, Pi], R2 = [Pi, p 2 ],---,R% = [Ps, Pmax], we prepare {3, 3, 3, 3, 3, 3, 3, 8} grid 
points. The distribution of the mesh points and the integration over each region are the 
same as those for the Amado model. 

It is a rather complex numerical task to get accurate solutions of Eq.flHSJ). We check our 
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numerical accuracy by reproducing the following optical theorem within 1% 
An 

— Im[t M B,MB(0 = 0)] = VMB,M'B> + CTMB,TmN (102) 

fc M'B'=itN,r)N,~fN 

where MB = ttN, rjN, yN are stable particle channels, the cross sections a a ,b are calculated 
from the non-resonant amplitudes tMB,M'B' by solving Eq. (}55|) . The two-pion production 
cross sections ctmb^-kN are calculated from the amplitudes Eqs.pH |) - pi|) with resonant am- 
plitude t^ BM , B , = 0. 



VI. RESULTS 

Our main interest in this paper is to use the numerical methods described in section 
V to examine the dynamical consequences of the one-particle-exchange interaction Z^ A7rA , 

ZnN,wA> an d ^-jv>A (FigJHI) • As illustrated in Fig|T3]and discussed in section IV, the ma- 
trix elements of these interactions have logarithmically divergent structure due to the 7T7tN 
unitarity cuts which are not accounted for in all of the recent calculations of two-pion pro- 
duction. The parameters needed to evaluate the partial- wave matrix elements of Z^ nA , 

ZpN,KAi an d ^v-H-A are fixed by the fitting the low-energy ixN and hit scattering partial- 
wave amplitudes, as given in Appendices D and E. With the resonant amplitudes also fixed 
by using the information of PDG to evaluate Eqs. ffUU|) - ff72jl . our first task is to choose the pa- 
rameters of starting Lagrangians, given in Appendix A, to evaluate the partial-wave matrix 
elements v[T s , M i B i LSMB {k' ^k) defined in Appendix C and v[T s , M , B , t \ \ N (k' , q) in Appendix 
F, with MB,M'B' = irN, r]N, 7rA, pN, aN. Here we are guided by the previous works on 
meson-exchange models of hN and NN interactions, as discussed in Appendix A. We also 
need to regularize the resulting matrix elements of all of the non-resonant interactions given 
explicitly in Appendices C and F. This is done by multiplying each strong interaction ver- 
tex in the considered non-resonant mechanisms, illustrated in FigsOH by a form factor 
[A 2 /(A 2 + k 2 )] 2 with k being the momentum associated with the meson at the MBB vertex 
or the meson being-exchanged. We adjust the cutoff parameters A as well as some of the less 
well determined coupling constants to get a reasonable description of the Jlab data of invari- 
ant mass distributions of yp — > 7r + 7r~p reactions. With the parameters listed in Tables I-II 
of Appendix A, our results (solid curves) of the invariant mass distributions are compared 
with the data at W = 1.88 GeV in FigUTA While the improvements are clearly needed, 
the chosen parameters are sufficient for our present very limited purposes of investigating 
the effects due to 7T7tV cut. No attempt is made here to adjust the parameters to fit all of 
the available data of yp — > tt + tt~p, tt°tt°p, 7r + 7r°n. This can be meaningfully pursued in a 
coupled-channel approach only when the data of ttN — > ttN, r]N, ttttN and 7V — > ttN, r)N 
are also considered. Here we focus on the effects due to the 7nrN cut which are neglected in 
all recent two-pion production calculations.. 

To see the dynamical content of our calculations, we also show in Fig ^] the contributions 
from each of the unstable 7rA, pN, aN channels. The M n + p distribution (top panel) is clearly 
dominated by the process •yp — > nA — > nuN (dashed). The peak near M w + p ~ 1.23 GeV is 
dominated by the ■yp — > it~(A ++ — > n + p) process, while the shoulder in the M w + p ~ 1.4 — 1.6 
region is due to the yp — > 7r + (A° — > n~p) process. The contributions from the pN (dotted 
curve) and aN (dot-dashed curve) are sizable and can change the shape and magnitude of 
the cross sections through interference effects. The M w + W - distribution (middle panel) is 
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dominated by ■yp — > p(p — > 7r + 7r~) (dotted) and hence is peaked at M n + n - ~ 0.76 GeV. 
However the contribution from ttA channel (dashed) are clearly important in getting the 
good description of the data. The situation for M^-p distribution (right) is similar to that 
for M n + p distribution (bottom panel), except that the relative strength between two A peaks 
is changed. 

We now turn to investigating the effects due to the one-particle-exchange driving terms 
Z 7rA ^ 7TA (E),Z pN ^ A (E), and Z^ vA (E) which contain the effects due to the nnN unitarity 
cut, as discussed in section IV. Their singularity structure, illustrated in Fig|T2*l is similar 
to that shown in FigEJof the three-boson case. We thus expect that non- resonant partial- 
wave amplitudes associated with ttA, pN, and aN states have similar momentum-dependent 
structure of FiglTol This is confirmed in our calculations. Some of our typical results are 
shown in FigJT7| for t nA>7rN and FigJTH] for the photo-production amplitudes t vArjN (upper 
panel) and t P N n N (lower panel). The solid curves in these figures are from our full calcula- 

( F) ( F) 

tions, which show rapid varying structure. When the driving terms Z^ A ^ A (E), Z pN ^ A (E), 

and Z^J A (E) are turned off in solving Eq. (jo3|) . we obtain slow varying dashed curves. 
Here we note that the momentum variable k in Figs JTTllTHl is related to the sub-energy 
a(k, E) = E — E s (k) for the resonant particle (A or p) to decay in the presence of a spec- 
tator particle s (n or N) with energy E s (k). Thus the full curves in Figs ll7tfTSl also reflect 
the rapid dependence on the sub-energy a(k,E). We emphasize that the rapid dependence 
of these amplitudes on the sub-energy a(k, E) is a necessary consequence of tttiN unitarity 
condition, as discussed by Aaron and Amadoj^], and is similar to what can be seen in the 
iiNN studies0, 0, [ti[ - Our results clearly indicate that the usual tree-diagram approx- 
imation should be used with cautions in interpreting the extracted N* parameters. The 
rapidly varying structure associated with an unstable particle channels must be taken into 
account in any phenomenological extraction of the partial-wave amplitudes. These were not 
taken into account in the early partial- wave analyses fff^ of the data of ttN — > mtN . 

If we further turn off the multiple scattering mechanisms in solving coupled- channel 
equation Eq. (jo3|) . we get the dot-dashed curves in Figs lT7IITHl The large differences between 
the dash-dotted curves and the solid curves indicate the difference between the dynamical 
coupled-channel approaches and the recent tree-diagram models. 

( F) ( F) 

We next examine the effects of the one-particle-exchange terms Z^ AnA {E), Z pN ^ A (E), 

and Z^ nA {E) on the differential cross sections of yp — > tt + tt~p. Here we set p as the outgo- 
ing 7r + momentum, q the relative momentum between 7T~ and p. Two of our typical results 
of the dependence of the differential cross sections da/(dM n - p dQ p dQ q ) on the azimuthal 
angle <fi of q are shown in Fig^Jwith the final it + tt~p kinematics fixed at M n - p = 1.23GeV, 
cos^p = 0.183, 4> p = —3.1 rad., and cos^g = 0.80(/e/t), 0.1S3(right). Our full results are 

the solid curves. The dotted curves are obtained when Z^ nA , Z p ^ nA (E), and Z^]^ A (E) 
are turned off in solving the coupled-channel equation Eq. (jo3j) . Clearly, the effects due to 
these one-particle-exchange terms are very pronounced in changing both the shapes and 
magnitudes of the differential cross sections. Similar results are also seen in our calculations 
for other values of p of the outgoing 7r + and q of the relative momentum of the outgoing 
7i~p system. The results shown in FiglT^l further indicate that the rapid varying structure 
of the amplitudes shown in Figs fTTTfTHI must be accounted for in any analysis of two-pion 
production. 

In the recent studies of two-pion production, the data of invariant mass distributions 
da/dM n N and da/dM n7T of yN — > ttttN are most commonly used to extract N* parameters. 



29 






M _ ( GeV) 



FIG. 16: The differential cross sections of "yp — ► ir + Tr p reaction with respect to the invariant mass 
M^+p (top), Mk+k- (middle), and M n - p (bottom) at W=1.880 GeV. The data are from Ref.0. 
The solid curves are from full calculations, The contributions from it A (dashed), pN (dotted) and 
aN (dot-dashed) to the invariant mass distributions of "fp — > it + it~p are also shown. 
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FIG. 17: The half-off-shell amplitudes nN(k, ko, E). The invariant mass of the outgoing A is 
1.232 GeV and the total energy is E=1.880 GeV. The left (right) hand sides are the real (imaginary) 
parts of the amplitudes with ttN in Sn (top) and P31 (bottom). The partial- wave quantum 
numbers for ttA state are indicated in each figure. The solid curves are from full coupled-channel 
calculations. The dashed curves are from setting Z^ AitA {E) = Z^ nA {E) = Z^ wA (E) = 0. The 
dot-dashed curves are from further setting multiple scattering terms of Eq.(55) to zero; i.e. setting 

tlTs'M'B' lsmb(^' ' 1 &i ^0 = v ijS'M'B' lsmb(^' ' 1 ^> Note that the matrix elements of v w a,wN are 
real in our phase convention (see Appendix A) and hence there is no dot-dashed curves in the right 
hand side. 



Since these cross sections involve integrations over angles of outgoing particles, as seen 
in Eq.(|73)l. the rapid varying structure of the partial- wave amplitudes due to nnN cut is 
washed out. We thus see the smooth distributions da/dM^N and da/dM nn , as shown in 
Figs|H3 However the one-particle-exchange terms Z^ A \ A (E), Zp^ wA (E), and Z^j %A (E) 
can change their magnitudes and shapes significantly. One example is shown in in Figl20l 
for 7£> — > -K°TT°p. We see that when these one-particle-exchange driving terms are turned off 
in solving coupled-channel equation Eq. (J55|) . the predicted invariant mass distributions are 
reduced significantly. Such a large difference further indicate the importance of including 
the mrN cut effects in calculating these particle-exchange terms for analyzing the two-pion 
production data. 



VII. SUMMARY AND FUTURE DEVELOPMENTS 



For analyzing the meson production data in the nucleon resonance (N*) region, we have 
developed a dynamical coupled-channel reaction model. With the assumption that the ba- 
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FIG. 18: The half-off-shell amplitudes t n & t yN(k, q, E) (upper) and t p N,-yN(k, q, E) (lower). The 
invariant mass of the outgoing A (p) is 1.232 GeV (0.76 GeV) and the total energy is E=1.880 
GeV. The partial-wave quantum numbers for the final irA and pN states are indicated in each 
figure. The solid curves are from full coupled-channel calculations. The dashed curves are 
from setting Z^ wA (E) = Z^^^{E) = Z^J vA (E) = 0. The dot-dashed curves are from fur- 
ther setting multiple scattering terms of Eq.(65) to zero; i.e. setting f^ MB A 7 Aivm r {k',k,E) = 
v lsmb x \ N m T {k',k,E). Note that the matrix elements of v^a^n { v pN^N ) are pure imaginary 
(real) in our phase convention (see Appendix A) and hence there is no dot-dashed curves in the 
right (left) sides of the upper (lower) parts. 



sic degrees of freedom of the considered reactions are mesons (M) and baryons (B), our 
starting point is an energy-independent effective Hamiltonian which is derived from a set of 
Lagrangians by using a unitary transformation method. Within the constructed Hamilto- 
nian, the N* excitations are defined by bare N* — > MB, ttttN vertex interactions and the 
non-resonant meson-baryon interactions are defined by the tree-diagrams generated from the 
considered Lagrangians. We then apply the standard projection operator techniques |6J] to 
derive coupled-channel equations for calculating the amplitudes of meson-baryon reactions. 
The model satisfies the unitary conditions within the channel space spanned by the con- 
sidered two-particle meson-baryon states and the three-particle nirN state. In this paper, 
we present explicit formulations within a Fock-space spanned by the basis states 7 A", tcN, 
t]N, 7T A, pN, aN, and tttcN. However, the formulation can be straightforwardly extended 
to include other meson-baryon states such as Kaon-Hyperon (KY) and uN, and other two 
meson production channels such as rjir N and KKN . 

To facilitate the interpretations of the extracted A^* parameters, we cast the reaction 
amplitudes into a form such that the meson-baryon scattering effects on A^* excitations 
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FIG. 19: Differential cross sections of 7p — * 7r + 7r~p at the jN invariant mass W=1.880 GeV. The 
outgoing 7r + momentum is p and the relative momentum between 7r~ and p is q. (j) is the azimuthal 
angle of q. The results are for the invariant mass M^- p = 1.23 GeV, cos 9 P = 0.183, <p p = —3.1 rad. 
The left (right) panel is for cosd q = 0.80(0.183). The dashed curves are obtained when Z m'b 
term is turned off in solving Eq.(65). 




FIG. 20: The invariant mass distributions of "yp — ► tt tt°p reaction at W=1.88 GeV. The dashed 
curves are obtained when Z MB M , B term is turned off in solving Eq.(65). 

can be explicitly calculated. These effects, called the meson cloud effects, are due to the 
mechanisms that the incident meson interacts with the baryons through all possible non- 
resonant scattering before the iV* is excited by the bare iV* — > MB vertex interaction 
of the model Hamiltonian. The determination of the meson cloud effects from the meson 
production data could be useful for interpreting the extracted iV* parameters in terms of 
hadron structure calculations. For example, it was found in Refs. (loL Hll ] that the meson 
cloud effects can account for the main differences between the extracted 7iV — > A (1232) 
resonance transition form factors and the constituent quark model predictions. It will be 
interesting to explore how the meson cloud effects, as defined in our formulation, can be 
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related to the current Lattice QCD calculations. 

In addition to giving a complete presentation of our theoretical framework, we also present 
in this paper a numerical method based on a spline-function expansion for solving the 
resulting coupled-channel equations which contain logarithmically divergent one-particle- 
exchange driving terms. These driving terms contain the effects due to the irnN unitarity 
cuts which must be included accurately in calculating the two-pion production observables. 
We explain how this method can be applied in practice for a simple three-boson Amado 
model, and then for our realistic model with 7 AT, tcN, ttA, pN, aN, and ttttN channels. 

An another important step in carrying out numerical calculations is to find an efficient 
way to calculate a large number of partial-wave matrix elements of the considered non- 
resonant meson-baryon interacting terms which are needed for solving the coupled-channel 
equations. Here we make use of the helicity representation of Jacob and Wick and also 
introduce a helicity-LSJ mixed-representation which is most convenient for calculating the 
electromagnetic matrix elements. While these are rather technical details, but are also 
presented explicitly in this paper for the completeness in explaining our numerics. 

With the parameters of the model chosen appropriately to fit JLab's two-pion photo- 
production data, we apply the developed numerical methods to show that the logarithmically 
divergent one-particle-exchange driving terms in the constructed coupled-channel equations 
generate rapid varying structure in the matrix elements of reaction amplitudes associated 
with unstable particle channels ttA, pN, and aN. Our results confirm the analysis by 
Aaron and Amado 0. We further show that these one-particle-exchange terms have large 
effects in determining the two-pion production differential cross sections both in shapes and 
magnitudes. Our findings suggest that one needs to be cautious in interpreting the iV* 
parameters extracted from the approaches which do not account for the effects due to the 
7nrN unitarity cuts. 

The calculations presented in this paper are far from complete within our formulation, 
while they are sufficient for testing the accuracy of our numerical methods and illustrating 
the importance of irnN unitarity cut. The N* parameters can be convincingly extracted 
and properly interpreted only when we apply our full formulation to analyze all available 
data of meson production reactions. Obviously this is a rather complex process. We now 
discuss how we will accomplish this rather ambitious research project in practice. 

Our first task is to fit the ttN elastic scattering data to fix the parameters defining the 
strong interaction parts of the model Hamiltonian Eqs.(j7j)- (|lUj) . This must be done by 
extending the coupled-channel calculations described in section IV in two aspects, First, we 
must include the driving term 4/m'b 1 defined by Eq.([3*T|). As discussed in section III, 
this term contains the irnN cut effects originated from the nNN vertex. Our second main 
task is to develop appropriate parameterizations of the bare iV* — > MB form factors for 
calculating the resonant amplitudes rigorously according to Eqs. (fTo'j) - (f*rT?|) . Here we need to 
make use of the predictions from hadron structure calculations. For example, we at least 
can fix the relative phases between different N* — > MB transitions by using the naive SU(6) 
quark model with meson-quark coupling. Predictions from more sophisticated models, such 
as the 3 Po model of Ref. [T^l and the model based on Dyson-Schwinger Equation could 
provide useful information to our investigation. In fitting the irN elastic scattering data, we 
should also fit the available ttN —>■ rjN reaction data and use the optical theorem to make 
sure that the predicted ttN total cross sections are also in agreement with the data. 

Once the irN data are fitted by the above procedures, most of the strong interaction 
vertexes in the non-resonant electromagnetic interactions f 7 iv,jwB an d v^n^tt-kN of our model 
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Hamiltonian have also been determined. We thus can focus on the determination of 7iV — > 
N* form factors. From Eq. (jl9|) . one can use the operator relations Eqs.(B32)-(B33) of 
Appendix B to write the dressed N* — > ^yN vertex of the resonant amplitude fEq. (jl5|0 as 



(E) — rjv*-» 7 iV + J! ^ N*-+MBGMB(E)tMB,iN(E) 



MB 



= ^N*^-yN + X! Fn*->MbGmb(E)vmB,"/N- (103) 
MB 

Since Tn*_mb m the second line of the above equation has been determined in the fit to 
the ttN reaction data, the bare N* — > ^yN vertex rjy.^ 7 Ar is the main unknown and can 
be determined by fitting the data of photo-production and electro-production of 7T, rj and 
two pions. Of course some less well-determined parameters in the non-resonant interaction 
v-yN,M'B' should also be adjusted in the fits. In practice, one can extract bare N* — > 
form factor at each Q 2 . It of course will be more interesting if the parameterization of the 
bare form factor Tn*^n can be guided by some theoretical calculations. 

We now turn to discussing the extension of the model to include KY and uoN channels 
which are also useful in probing the structure of N*. In particular, we note that r]N, KA, and 
uoN channels are of isospin T = 1/2. The properties of T = 1/2 N* states can therefore be 
more selectively extracted from analyzing the production data of these three channels. Thus, 
an extension of the formulation presented in this paper to include KY and ujN channels 
is highly desirable and technically straightforward. However, it will increase the needed 
computation effort enormously. Nevertheless, we can make use of the results from fitting 
the irN, rjN and ttttN data to perform simplified coupled-channel analyses of the KY and 



ujN production data. This can be done by following the approach of Ref. |81 . 

Considering the KY production, we assume that it can be described by a coupled-channel 
model including 7 AT, KY, ttN, and a dummy channel QQ which represent all of the neglected 
channels. If we further assume that KY does not couple directly with the QQ state (mainly 
because there is no information about how KY couples with ttttN channels), one can cast 
the coupled-channel equation Eq. ()24|) into the following form 



with 



t~iN,KY(E) = v iN)K y\\ + G KY (E)t KY ,KY(E)] + v yN ^ N G 7rN (E)t 7rNjKY (E) (104) 

tKvM E ) = v e J YiKY (E)[l + G KY (E)t KY , KY (E)}, (105) 
t KY M E ) = \l + t KYtKY {E)G KY {E)}v KY ^ N {l + G^ N {E% N ^ N {E)}. (106) 

Here the effective KY interaction is defined by 

v ky,ky{E) = v KYjKY + v KYnvN G„ N {E)[l + U Ni7 , n (E)G wN (E)]v wN! ky , (107) 

and t-jxN^N is from solving the coupled-channel equation Eq.(j2U) in the ttN © QQ space. 

If we assume that the dummy channel QQ = r]NQ)iTAQ)pNQ)<TN, the scattering amplitude 
tirN^N m the above equations is just the solution of Eq.(j21J) of the model determined in the 
fit to 71 N data described above. We therefore can use this information to solve Eqs. (|105|) - 
()107|) and determine the parameters associated with the non-resonant interaction vky,-kN 
and vky,ky by fitting the available data of nN — > KY reactions. This will then allow us 
to generate txY,KY and t n N t KY to evaluate Eq. (jl04|) and also fix the strong vertexes in the 
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non- resonant v^n^y- The KA photo-production and electro- production data can then be 
used to extract the 'yN — > N* form factors for T = 1/2 iV* states. The same procedure can 
be used to analyze the uN production data. 

To end this paper, we would like to emphasize here that the objective of performing 
dynamical coupled-channel analyses of meson production data is not only to extract the 
N* parameters, but also to provide information on reaction mechanisms for interpreting the 
extracted N* parameters in terms of the quark-gluon substructure of hadrons. In particular, 
we account for the dynamical consequences of the ttttN unitarity condition which is very 
difficult, if not impossible, to be treated rigorously in the existing approaches for calculating 
the hadron structure or the Lattice QCD calculations. An another important point to note 
is that our approach accounts for the off-shell scattering effects which describe the meson- 
baryon scattering wavefunctions in the short range region where we want to explore the 
structure of N*. These essential quantum-mechanical effects are absorbed in the parameters 
of the approaches based on tree-diagram models or K-matrix models. Thus our dynamical 
approach perhaps has a better chance than these two approaches in revealing the quark-gluon 
substructure of baryons. Our progress in this direction will be published[67] elsewhere. 

We would like to thank B. Julia-Diaz and K. Tsushima for their help in checking our 
calculations of the matrix elements of non-resonant interactions. This work is supported 
by the U.S. Department of Energy, Nuclear Physics Division, under Contract No. W-31- 
109-ENG-38 and the Japan Society for the Promotion of Science Grant-in- Aid for Scientific 
Research (C) 15540275. 
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APPENDIX A: LAGRANGIAN 



In this appendix, we specify a set of Lagrangians for deriving the non-resonant interactions 
vmb,m'b' which is the input to the coupled-channel equations Eq.pijl. Here we are guided 
by the previous works on meson-exchange models of ttN and NN interactions. The coupling 
with pseudo-scalar mesons 7r and 77 are consistent with chiral symmetry. The vector meson 
couplings are less known and are mainly constructed phenomenologically. In the convention 
of Bjoken and Drell|82|), the Lagrangian with 7r, 77, N, and A fields are 

' / l ^ N -i>NlvT{hTil) N ■ , (Al) 



j ttNA 



h VA -^f ^ • , (A2) 



L^AA = -^-fpA^^TA^A ■ 9u$tt , (A3) 
L V NN = -^ 1 ^-i ) Nlnlb^Nd IJ, (j)r 1 - (A4) 



The interactions involving p meson are 



LpNN = QpNN^N [^ffM ~ 7r £ — a ^ dU ]fi 1 ' n^N , (A5) 

LpNA = -i^^ i ipAl u l5T ■ [dppl - dvp^tpN + [/i.e.] , (A6) 

Tflp 

LpAA = gpAA^Aal^ ~ ^^d u ]p^ ■ , (A7) 

Lpw = gpmlfar X dpfa] (A8) 

LnN P tv = gpNN^NlyTi^N ' X </>tt , (A9) 

2 

L NNp p = -^^^n^tiPn ■ p-p x pi . (A10) 
8m N 

Note that the contact terms Eqs. (jA9j) - (jA10|) are from applying [<9 M — > <9 M — gpNNp 1 *-} on 
L nNN Eq.flM) and L pNN Eq.flMJ. 

The interactions involving u meson are 

Lunn = 9uNNipN[lp - —^-cr^d^oj^jsi , (All) 

2m N 

L^p = - 9 -^ep aXu d a p»d x ^ . (A12) 
We also consider interaction involving a scalar isoscalar a meson 

LaNN = gaNNi>N^N4>a (A13) 
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0.08 2.049 1.825 11.5 38.4329 
TABLE I: Coupling constants determined in Ref . [lo| . 

g P iriT g<TiriT gmnp 9pNN fr/NN QcrNN firAA fpNA fpAA K pAA 

6.1994 1.77 11.2 6.1994 1.77 12.8 1.78 -6.08 -4.30 6.1 
TABLE II: Coupling constants used in the calculations in this paper. 

£»w = -f^^MAA- (A14) 

To proceed, we need to know the coupling constants of the above Lagrangians. The 
parameters determined from fitting the irN data within the SL model [10( are given in Table 
I. The p,a —>■ 7T7T coupling constants can be estimated from fitting tttc phase shifts in the 
isobar model[68j, as described in Appendix C. The decay width of uj — > up can be used to 
estimate the coupling constant g^ p . The rjNN coupling constant f V NN has been estimated 
in recent studies of i] production from irN and 7 A" reactions. The aNN coupling can be 
estimated from the previous works on AW scattering. These parameters are adjusted around 
the values from these estimates to fit the JLab data of 7]? — > tt + tt~p reactions, as described 
in section VII. They are listed in Table II. 

We have very little information on the coupling constants f^AA, Ina p and f p AA- We 
simply follow the previous works and use the simple SU(6) quark model to determine them 
from the empirical values of the coupling constants f^NN and q p nn- To be more informative, 
we here also describe how this procedure is used in practice. 

First step is take the static-baryon limit of the matrix elements < B' \Lmbb 1 \BM(q) > 
to define the effective MBB 1 Hamiltonian operators in the spin-isospin space of baryons. 
They are 

H^nn = i^^-a ■ qr a , (A15) 

m n 

H^na = i^S-qT a , (A16) 

HnAA = z^^A-gT^, (A17) 
m n 3 

HpNN = i g o NN } 1 + ^ Z x q ■ e(p)r« , (A18) 

4:171 N 

HpNA = -i fj ^S x q ■ e(p)T a , (A19) 
HpAA = -igpAA 1 + 2 ^ pAA ls A x q- € p Tl . (A20) 

Here, a is the isospin component of the considered meson, 5* and T are the spin and isospin 
operators of the N-A transition, Sa and 7a are the spin and isospin operators of the A. 
Along with the usual Pauli operators a and r, they are defined by the following reduced 
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matrix elements 



< N\\a\\N > = < N\\t\\N >= y/6, (A21) 

< A||S||JV > = < A||T||JV >= 2, (A22) 
< A 1 1 S" A 1 1 A > = < A||T A ||A >= y/lE (A23) 



with the convention that 



< jfrnflO^ljiTm >= 1 < jfmfljilrriiM >< jfWO 1 ^ > . (A24) 



12 j, + 1 

We next consider a simple meson-quark interaction Hamiltonian 

H, m = f -^Y. & -fa, (A25) 

m i" i=l,3 

H pqq = ^ £ i*i x 9 • e > > ( A26 ) 

m P i=l,3 

where cr^Xj are the spin and isospin operators of the constituent quarks. By using the Os 
constituent quark wavefunctions ipN,m aN m TN and ^.m,,^ f° r the nucleon and A and the 
relations Eq. (jA21|) - (jA24|) . we have the following relations between the matrix elements in 
the spin-isospin space 

5 

< ^N,m' 3N m' TN \ °i T i\^N,m SN m TN > = ~ < m'^KJ^l^K,, >> ( A<27 ) 

i=l,3 6 

< i)A,m> SA m> TA \ £ < ^N,m SN m TN > = 2V2 < m' SA m' TA \ST\m' SN m! TN > , (A28) 

i=l,3 

4 

< ^A,m' 3A m' T J £ ^ r #A,m SA m TA > = o < m s A m t A I 5aTa l m *A m T A > . (A29) 

i=l,3 

Using the above formula and assume that the matrix elements of the hadron Hamiltoni- 
ans Eqs. (jA15)) - (jA20|) are equal to the matrix elements of the quark- meson Hamiltonian 
Eqs. f)A25|) - (jA26|) within the SU(6) chiral constituent quark model 

< ^B, m ' SB , m ' rB \H Mqq \ip BtmsBimTB >=< m' SB ,m' TB \H MBBI \m SB ,m TB > , (A30) 
we then obtain 



fnNN — 


If 

g J nqq ) 


(A31) 


firNA = 




(A32) 


-JttAA — 


if 

g J nqq i 


(A33) 


fpNN = 


If 

2 J pqq ) 


(A34) 


fpNA = 


2\/2 f pqq , 


(A35) 


g/pAA = — 


2 J pqq j 


(A36) 
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where we defined 



. 9 p nn(1 + k p) 

JpNN — : wn P , {-^3 1) 



fpAA = 9 P AA l+ n KpAA m p . (A38) 
zrfiA 

From the above relations we finally have 

[72 

fnNA = y —firNN , (A39) 

6 

fnAA = ^firNN , (A40) 
[72 

fpNA = -]J-^fpNN, (A4l) 
6 

fpAA = ^ fpNN ■ (A42) 

By using the vector meson dominance assumption and the recently determined A magnetic 
moment, we can set 

k p aa = 6.1 . (A43) 

With the values in Eq. (|A43|) and values listed in Table I, we can use Eqs. flA37|) - (jA42|) to get 
fnAA, fpNA and fpAA- The resulting values are also listed in Table II. 

The electromagnetic interactions are obtained from the usual non-interacting Lagrangian 
and the above interaction Lagrangian by using the minimum substitution dp — » dp — ieAp. 
The resulting Lagrangian are given below : 

L jN n = Me N Y - T^cr^d^NAp , (A44) 

V = [4 x d"^] 3 Ap , (A45) 

L^nkN = ^^-[^nYi^n) x (f^hAp , (A46) 

Lypp = Wp» - <9V) x p„} 3 Ap , (A47) 

L ypmt = -gpnAiP^ 1 x &r) x ^hAp, (A48) 

L iN ,a = — [(^^Var) x ^]sAp } (A49) 



L~/NpN 


— 9pNN[ 


L-yNA 


= -i^A- 


-LiypTr 


= e< 






= L e, 


m n 


L-ypr] 




m p 



2rriN 2 



(A50) 
(A51) 

( ^4-(dV)(^), (A52) 
^{d^A^lidnJ) , (A53) 
v ^d^ v d a A^ (A54) 
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= -^M(aM"-9M>, (A55) 

Tflp 

L-yAA = r A (T 3 A + l -)[-i»g w + « 7 * + rf-f) + Iiv^WaA, • (A56) 
For Eq. (| A44j) . we have defined 

„ F 15 + F iy r 3 



2 

F 2S + F 2V t 3 



(A57) 
(A58) 



where F 15 (0) = F 1V (0) = 1, F 2S (0) = ^ p + p n - 1 ~ -0.12 and F 2y (0) = /i p - fi n - 1 ~ 3.7 
The matrix element of 7AA vertex of Eq. (jA49|) between an N with momentum p and a A 
with momentum p& can be written explicitly as 



< a(p' a )\t^ a \n( Pn )> 



m A + m N 



^ 2m N (m A + rn N ) 2 — q 2 

12 

+G E il5 7 ^e llXa pP a q P e x uaS p\q S 

(m A — mAr)^ — q 2 

+GcijB-, ^{q 2 P v -q-Pq v )i (A59) 

(m A — m N ) z — q z 

with P = (p' A + Pn)/2 an d p' A = p^ + q. Note that the index p of contracts with the 
A field and v with the photon field. The coupling strength Gm = 1-85, G E = 0.025, and 



Gq = —0.238 are taken from the SL modelllOj, 111 



APPENDIX B: DERIVATION OF COUPLED-CHANNEL EQUATIONS 

In this appendix, we give the derivation of coupled-channel equations from the model 
Hamiltonian H e ff = H Q + V defined by Eqs.(|5|)-([TU|). We apply the standard projection 
operator techniques . The procedure is similar to that used in the derivation of nNN 
equations (65l|. We start with Eq.(jTJJ 

T(E) = V + V-^— T(E) . (Bl) 
£/ — li 

The propagator in the above equation is understood to include +ie for defining the boundary 
condition, but is omitted to simplify the presentation in this appendix. The interaction V, 
defined in Eqs.(l7j)- (frUj) . can be more clearly written as 

V = V 22 + V 33 + (T 12 + U 23 + 713) + (r 2 l + ^32 + 73l) , (B2) 

where v 22 = v M b,m'b> + v n7r , T 12 = T N ,^ MB + h M *-^n with M* = p,a, 713 = T N *^ n7T N, 
u 23 = Vmb,-kttNi and v 33 = v^Nt^n- Here we restrict MB = jN,7rN,r}N,7rA, pN,crN. In 
Eq. (jB2|) . we have also introduced more transparent notations T 21 = r{ 2 , u 23 = u 23 , and 

731 = 7l3- 
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We next introduce projection operators 

P + Q = 1 (B3) 



with 



Q = 1 7i7i N >< 7i7iN\ , (B4) 
P = P x + P 2 + P 2 * , (B5) 



where 



P 1 = J2\ N * ><N*\, (B6) 

TV* 

P 2 = | 7 jV >< 7 iV| + |vriV >< vriV| + \r]N >< r}N\ , (B7) 

P 2 „ = |ttA >< ttA| + \pN >< pN\ + \aN >< aN\ . (B8) 

We then obtain the equations for the projected operators T PP = PTP and Tq P = QTP 

Tp P = Vpp + Vpp— tj-Ppp i (B9) 

— li 

T QP = -i t7 ~ q -VqpI 1 + -^V^pp] , (BIO) 
1 " QQ e-Hq-Vqq E ~ H » 



where 



with 



Vpp = Vpp + V PQ g V QP , (Bll) 

£j — tlQ— Vqq 



Vpp = pvp = v 22 + r 12 + r 21 , (Bi2) 

V QP = QVP = u 32 + r 21 + 73i , (B13) 

Vqq = QVQ = v 22 + v 33 . (B14) 

Eq. ()Bll|) can be written explicitly as 

Vpp = p[(v 22 + r 12 + r 21 ) + (u 23 + 7l3 + r 12 )G Q (M32 + 731 + r 21 )]p , (bis) 

with 

GQ = E-H -Q(v 22 + v 33 )Q ' (B16) 

From the definitions Eqs. (jB6| )- (jB8|) for the projection operators, we have the following con- 
ditions 

P 2 r 12 Q = gr 12 p 2 = o, 

P2731Q = P 2 *73iQ = Q131P2 = QlziP2* = . (B17) 
With the above "doorway" conditions, we can decompose Vpp as 

Vpp = P[S + v)P , (B18) 
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where 

£ = [Ti 2 GqT 

2l] un— connected ■ 

(B19) 

Here un — connected in Eq. (jB19|) means that the pion emitted from one baryon is also 
absorbed by the same baryon. Obviously this is the self-energy of the unstable particles in 
the Ti A, pN and aN states of P 2 * space. We thus have 

PEP = P 2 *EP 2 * . (B20) 

All other interactions within the P-space are in v of Eq. (|B18|) 

v = V E + f 12 + f 21 + £ (B21) 

with 

V E = v 22 + (u 23 + T 12 )G Q (u 32 + T 21 ) - E , (B22) 

r 21 = r 21 + +Ar 21 , (B23) 

Tl2 = Tl2 + +Ar 12 , (B24) 
S = 7i3G Q 73i , (B25) 

where AT 2 i and Af \ 2 contain interactions due to iV* *-> irirN transitions 

Af 21 = [u 23 + r 21 ]G Q 73i , (B26) 
Af 12 = 713^^32 + T 12 ] . (B27) 

To follow the derivations given below, we note that the well known operator relations 

= v + t-±—v (B28) 
t, — Hq 

lead to 

= v[l - - 1 TT v}' 1 . (B29) 



E-H 



Eas. flB28|) and (lB29jl then lead to 



[i-v — - — r 1 = i + t- 1 



Ea. (|B28|) also leads to 



E — Ha E — Hq 



t = v + v E^T- v - (B31) 
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Comparing Eqs. (jB28|) and ()B31|) . we have 

1 1 



-t. 



E-Hq-v E-H 

"^db^'^iv (B32) 



It can also easily be seen that 



1 -t-Asr. (B33) 



E-Hq-v E-H E-H E-H 



In the following derivations, the above relations Eqs. ()B28|) - (jB33|) will be often used without 
mentioned them again. 

By using Eqs.(jB28}, (IB3T1) and (lB33j) . we can write T PP defined by Eq.flBH) as 

T PP = P[(j: + v ) + (j: + v ) 1 (E + v)]P 

hi — SlQ — Zj + V 

= P[S + S E-^o-S S + (1 + S E-ffo-S m(1 + E- P h]-^ )]P (B34) 

with 

T * = v + v E _^_ Y T,. (B35) 

By using Eq.(B13) and relation (B32), we can write Eq. (jB10|) as 

T QP = Q[(l + t Q -^—)(u 32 + 73i + r 21 )[l + -L—T PP }P , (B36) 
hi — riQ hi — iiQ 

where 

tQ = V QQ + Vqq — ZTff tQ ( B37 ) 

describes ttttN — > ttttN scattering through V QQ = Q[v 2 2 + v 33 }Q = v n7T + v nNj7rN + v n7TN ^ nN 
interactions. 

We now derive equations for calculating the scattering amplitudes between two particle 
channels in P' 2 = P 2 + P 2 * space. We first note that 

(B38) 
(B39) 
(B40) 
(B41) 

The above relations and Eq.(B35) lead to 

P%P' 2 = V E + V E Pp,!* + f 21 Pl Q P{T,P' 2 , (B42) 

hi — hlQ — Zj hi — T7l)y* 

PiTvPo = f 12 + f 12 ^ PoTvP; + t I -^P l TvP; . (B43) 

1 v 2 12 12 E-H -J: 2 v 2 E-m%* 1 v 2 V ; 



P 2 VP 2 


= V E 


P' 2 vP x 


= f 21 


P1VP 2 


= f 12 


PivPi 


= t. 
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Eq. (|B43|) can be written as 



1 P' 

PiT^p' = fi-s ^l _ Mri 2 + r 12 2 - PoTvP'} 

1 v 2 [ E-m%. 1 P-P -£ 2 2J 

E - m%. ^ , . P^ 



[r 12 + r 12 — — j — -p' 2 t- v p' 2 \ . (B44) 



Substituting Ea. (lB44l into Ea. (lB42l) . we have 



where 



P^ = X + X— -1 -P^ , (B45) 



X = V ' + t * E-W.-t F >'- <B46) 
Eq. (|B45J) can be written as 

P{T,P 2 = [1 - 1^ E _^_ E ]-'[V E 

= iE + ^-a^.-^-^U-^ " *° " E1 ' (B47) 



where 



t B = y E + V F - -t B , (B48) 

r 21 = [l + ** g _g o _ E ]r a i • (B49) 

We further note that 

1 1 



P - P - £ - X ^_iy _E-VE-r 



12 B-M°.-f: r21 



+ 5 F7 1 ^ F F7 1 ^ t / (B50) 



P - P - S - V E E - H - £ - Vb P - H - E - V, 



with 



*. = f 12 — — ^ g-f 21 [l+ - - 1 - — t s 



P-M°, -£ P-Pq-E-Vz 



12 



P - Mjy* - E - E 
45 



f 2i , (B51) 



where 



S " Tl2 E-H -Z-V E T21 

= ri2 jE _^_ s r 21 . (B52) 



Here r 21 has been defined in Eq. (|B49|) . 

By using Eqs. (|B50|) and ()B51|) . Eq. ([B47[) can be written as 

P' 2 T^P' 2 = t E + f 21 - — ^-[1 + (f 21 - ——-—f 12) 



E — Afjy, — £ v E - H - E - Vb £ - M° . - £ - £ J 
+ ^21— — ^ — ^12— — — [E - H - £] 



E - M%* - £ - E ^E-Ho-^-Vi 



E 



fE + r " g -Ma.-E-£ f " 11 + ¥3faM. 



The above then gives 



where (also recalling Eq. (|B49jl ) 

f a = f 12 [l + e _h ^ b] , (B54) 
T 21 = [l + ^ ^_^_ s ]f 21 . (B55) 

We now turn to deriving equations for calculating two-pion production. For initial ttN 
or 7iV of P 2 -space, Eq. flB36|) can be written explicitly as 

Tqp 2 = Q(l + t Q ® )[u 32 + u 32 P2 T PP P 2 

hi 11Q hi 11Q 

+ (u 32 + T 2l )P 2 *-J^—T PP P 2 + 73i F Pl u T PP P 2 ] . (B56) 
hi — hiQ hi — Hq 

From definition Eq. (jB34|) . we have 

P{r PP P 2 = P 2 T- V P 2 , (B57) 
P 1 T PP P 2 = P{T- V P 2 , (B58) 

P2*TppP 2 = P 2 *(l + E— — —P 2 *]T ij P 2 

hi — Hq — 2j 

= P 2 ,[E - H }- ]- -P 2 *T,P 2 . (B59) 

hi — Hq — 2j 
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By using the above relations and Eq. (jB53|) and P' 2 = P 2 + P2*, the 3rd in the bracket of 
Eq. (jB56|) can be written as 

3rd = (u 32 + r 21 )^§-^ + ^ E _ Mi ]_ ± _^ ■ (B60) 

By using Eq.(B45), the 4th term in the bracket of Eq. ()B56jl becomes 

*th = 731^^ E f 12 [1 + g _^ — f»r,fa 

731 - + Fib"' + ^E-A^-E-E^ 1 





-M%* 


- E 




1 








- E 




1 




E 


-M%* 


- £ 




1 





731 - - - lFl2(1 + fl-lo-S ^ + (fl2 E-^ -S r2l) E-M°!-S-S ri2] 

- 1 

731 1 + S s =r r i2 

1 , 

= 731 s -IV B61 

m E-M°, -E-E 1 ; 

We finally obtain 

Tqp 2 = Qn { ~Jl[u 2 3 + {U32 jp 2 U + (^32 + 1^21 )- ^ ^}{^ + f 2 l- — ^ ^f 12 } 

a* — -no t, — — 1^ E — Mft* — E — E 

+ ^-m^- S -s f " |a - (B62) 

where 

n ^ = ( 1 + *0^r^)- (B63) 

Here tq is defined by Eq. (jB37|) and hence fi^Ar * s the t 1 " 71 "^ scattering operator. 

In the above rather detailed derivations, Eqs.(B53) and (B62) are what we need to in- 
vestigate meson-baryon scattering and two-pion production. In practice, the interaction 
731 = r^v* — >7r7r n will be neglected in first calculations. If we set 731 = 0, we then find from 
Eqs.(B23)-(B27) that 

f 2i -> r 21 , (B64) 
fi 2 -> T 12 , (B65) 
E -> 0. (B66) 

Eqs.(B52), and (B54)-(B55) lead to 

f 2i - [l + ^ E _^_ s ]r 21 , (B67) 

r 12 - r 12 [i + E *s] , (B68) 

s ~" ri2 ^-i -s [1 + tB ^-io-s ]r21 

= V ^ E _ l H _ s F 2i. (B69) 
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Recalling Eq.(B7)-(B8) for the projection operators P 2 and P 2 *, we can write 

P! 2 = P 2 + P 2 * = E \ MB >< MB \ > ( B7 °) 

MB 

where MB = jN, 7rN, r)N, 7rA, pN, aN include all meson-baryon states in the considered 
model space. Defining 

T M b,m'B'(E) = < MB\P' 2 TvP' 2 \M'B' > , (B71) 
t M B,M'B'{E) = < MB\t E \M'B' > , (B72) 
V M b,m'b'{E) = < MB\V E \M'B' >, 

and 

Gmb(E) = < MB\— A — I MB > , (B73) 

using the simplifications Eqs.(B64)-(B69), Eq.(B22) for Ve, and Eq.(B48) for the matrix 
element of Eq.(B53) between two MB states then become 

Tmb,m'B'{E) = t M B,M'B'{E) +t MB)M , B ,(E) , (B74) 

where 

tMB,M'B'(E) = VmB,M'B'{E) + E VMB,M"B"{E)G M "B"{E)t M » B » tM i B i(E) , (B75) 

M"B" 

with 

Vmb,m'b' (E) =< MB\v 22 + (u 23 + T 12 )G Q (u 32 + T 21 ) - E\M'B' > . (B76) 

As defined in the beginning of this appendix, we have v 22 = vmb^'B'+v^, Ti 2 = Tn*->mb + 
h M *^ niT with M* = p, a, T 2 i = T\ 2 and u 23 = vmb^ttN, u 32 = u\ 3 - Eq. ()B76|) can be written 
explicitly as 

Vmb,m'B'(E) = v M b,m'b' + Z MB ^ M i B t(E) . (B77) 

Here Zmb,m'B'(E) contains the effects due to the coupling with titiN states. It has the 
following form 

Z M b,m'B'(E) = <MB\ F- — ^ —Ft I M'B' > 

E - H Q - v„tt N + xe 



where 



-[Smb,m'B'^mb(E)] , (B78) 



E MB (E) = < MB\E\MB > , (B79) 

F = 9V + VmB,ttttN 

= [r A — >7rN + h p ^ n7T + h/j — >777fJ + VMB,TTTtN , (B80) 

VmrN = V n N,nN + V wn + V n7r N,irnN ■ (B81) 



48 



The resonant term in Eq. (jB74|) is 

tuB,M'B>( E ) = <MB\T 21 ^—^—^Y 12 \MB> 



J2 Tmb^n* < N t\ E _l > Tn^m>b> ■ (B82) 



N* ,N* 



Note that E in Eqs.(B73) and (B79) is defined by Eq.(B19). If we neglect the non-resonant 
interactions in tttcN Q-space, we then have 

E ^ [1^12"^ 2l\un— connected • 

(B83) 

±j — n 

Since Ti2 does not have aiV^ ttN, we obviously have < tt^I S |7rA^ >= and hence 

° ME) = E-K„(k)-K N (j,) + ie' (B84) 

GUE) = E-KM-KM-X^E-K^k))' (B85) 

° ME) = E - K p {k) - K N (p) - E pN (E - K N (p)) ' (B86) 

GME) = E-K a -K N {p)-* aN {E-K N {p)y (B87) 



where 



1 

EttaM = <nA\T A ^ 7vN — — , T wN _> A \iiA > , (B88) 

uj - K 7T (k) - K N {p) + te 

Z pN (u) = < P N\r p ^ 1 V^ p \pN > , (B89) 

u - E^ki) - E 7T {k 2 ) + ie 

Y, uN {uj) = < aN\Y^ 1 , r^ a \aN > . (B90) 

a; - E n (ki) - E 7T (k 2 ) + ie 



In the above equations, K a (p) = \jm 2 a + p 2 is the free energy operator defined by momentum 
operator p. 

When N* — > tttiN is neglected, the two-pion production operator Tq P defined in 
Eq. (B62) also becomes simpler, since its last term in the right-hand side does not contribute. 
By using Eqs.(B66), (B72) and (B77), the matrix element of Eq.(B62) T^ N)MB (E) =< 
itttN\Tqp 2 \MB > can be written as 

T^ N , MB (E) = <ip { -J N \u 32 \MB > 

+ £ [< % \M'B' > Tm'B' ,MB 

M'B> & - n - 

+ < ^\ t ^ E - P hI-y} m ' b ' > Tm ' b '^ ■ ( B91 ) 

Recalling that u 32 = iWa^mb, r 2 i = T^n^a + r^-^ + r^^cr, we can write Eq. (jB91j) 
explicitly as 

TttkN,mb{E) = T^ NMB (E) + T%£ N)MB {E) + T^ NMB (E) + T^ NjMB {E) (B92) 
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with 



T% N>MB (E) = <4*n(E)\ E v™n,m'b>\M'B'>[6 m >b>, M b 

M'B' 

+Gm'B'(E)Tm'B',mb(E) , (B93) 

T:Xmb(E) = < ^-J N (E)\r nN ^\nA > G^{E)T^ MB {E) , (B94) 

T^ NMB {E) = < ^-J N (E)\K^ p \pN > G pN (E)T pN ,mb(E) , (B95) 

T^ N)MB {E) = < ^ N {E)\h^ a \aN > G aN (E)T aN>MB (E) . (B96) 

APPENDIX C: MATRIX ELEMENTS OF MESON-BARYON POTENTIALS 

To solve Eq. (|B3j) for generating the non-resonant amplitudes, we need to first calculate the 
partial- wave matrix elements of meson-baryon non- resonant interactions vmb,m'b' generated 
from the Lagrangians specified in Appendix A, and the one-particle-exchange interaction 
^mb M'B' (E) defined by Eq. (JHUj) and illustrated in Fig|SJ In this appendix, we present 
formula for calculating the partial- wave matrix elements of vmb,m'B' with MB, M'B' = 
irN, r]N, aN, pN, 7rA. The partial-wave matrix elements of Z^b m'B'(E) will be given in 
Appendix D. 

In general, each of the constructed vmb,m'b' consists of various combinations of tree- 
diagram mechanisms illustrated in Fig|21 They can be computed by the usual Feynman 
rules, except that the time components of the propagators of the intermediate states are 
specified by the unitarity transformation method, such that the resulting matrix elements 
are independent of the collision energy E of Eq.([55|) and free of any singularity on the real 
momentum axis. We will explain this feature of our model at the end of this appendix. 

It is convenient to get the partial matrix elements by first evaluating the matrix elements 
of v mb, M'B' i n helicity representation and then transforming them into the usual \ (LS)JT > 
representation with J, T, L, and S denoting the total angular momentum, isospin, orbital 
angular momentum, and spin quantum numbers, respectively. For each meson-baryon (MB) 
state, we use k(p) to denote the momentum of M(B). In the center of mass frame, we thus 
have p = —k. Following the Jacob- Wick formulation [83j] . the partial- wave matrix elements 
of the non-resonant interaction Vmb m 1 b 1 can be written as 



v 



JT 



; (k',k,E) = £ 



(2L + 1)(2L' + 1) 



L'S'M'B',LSMBV* i ^ ) — L 2J + I 

x < j'mJ'b^m - A' s I S'S' Z >< L'S'0S' z | JS' Z > 



x < JmJb\m - Xb I SS Z >< LS0S z \ JS Z > 

x < J, k'\' M - \' B | %'b',mb | J, k\ M - X B >) , (CI) 

where jm and js are the spins of the meson and baryon, respectively, and Xm and A# are 
their helicities, and 

< J, k'X' M — X' B | vm'B' mb I J, kXu — Xb > 
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Channels 


ttN 7]N 


aN pN 


7rA 


irN 


1 2 


4 7 


11 


rjN 


3 


5 8 


12 


aN 




6 9 


13 


pN 




10 


14 


7rA 






15 



TABLE III: Labels for v afi with a, b = ttN, tjN, aN, pN, and vrA 

= 2 7 r£ 1 rf(cos^)^- A ' fl , AM _ AB W 

X < (A/, s' M \' M ), (-k', S B , -X' B ) | V M 'B',MB I S A fA Af ), (— fc, s B , -A B ) > 

Here we have chosen the coordinates such that 

k' = (k' sm9,0,k' cos6) , 
k = (0,0, k) , 

and the helicity eigenstates are defined by 

k-s M \ M(k, s m Xm) > = \m \ M(k, s M X M ) > , 

[-k-s B ] | B(-k,s B X B ) > = \b \ B(-k,s B X B ) > ■ 



(C2) 



(C3) 
(C4) 



(C5) 
(C6) 



Note the — sign in Eq. ()("6j) . 

To evaluate the matrix elements in the right hand side of Eq. ()C2)) with the normalization 
defined by Eq.([48|l. we define (suppress the helicity and isospin indices) 



< k'(j),p' | v M 'B',mb I k(i),p > 



1 



B 



m B 



xu B/ (p')V(n)u B (p) (C7) 



where n defined the MB — > M'B' transitions as specified in Table ITTT] and i, j are the isospin 
indices of the mesons. We also have defined q = k! — k or q = p — p'. The expressions of 
each term in Table 11111 are given in the following subsections. 



1. 7r(fc,i) +N(p) — 7r(fc',j) +N(pf) 



V(l) = vl + H 1 + V} + + V} 



(C8) 
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with 



K 1 = [— ] 2 ^ 5 r j S N (p + k) ^ (C9) 

v. 1 = [^] 2 for'Ssip - k') for* (CIO) 

yi = \!iE±f ka (Ttysf( p -k'%T3 (Cll) 

1 ^ 1 

- WpNN9p7T7T—eju^ 2 

2 q z — mj 

x p+ ^) + t^{(^+ *o ^+ r)}] (ci2) 

4mjv 

K 1 = -9aNN^S iJ 4 1 ^2 ( C13 ) 



2. w(k,i) + N{p)^rj(k')+N{p') 

V(2) = V a 2 + H 2 (C14) 

with 

K 2 = AjVAf/??AfjV ^75^(P + fc) ^ (C15) 

V b 2 = f ^ NN ^ NN K l5 r*S N (p - k') ^75 (C16) 
m n m v 

3. v (k) + AT(p) -► r/(fe') + iV(p') 

t/(3) = K 3 + H 3 (C17) 

with 

K 3 = [^] 2 #V?;v(p + ^75 (C18) 
H 3 = [^] 2 ^iv(p - k>) jt' l5 (C19) 

4. 7r(M) + W(p) -»• <r{U) + 7V(p') 

1/(4) = K 4 + n 4 + K 4 (C20) 
52 



with 

K 4 = ig ff NN—S N (p + k) ^ 5 r l (C21) 



H 4 = ig*NN^ frt>S N (p - k'V (C22) 



f>4 .firNN9cmiT , j ? • fc /n n o\ 

K = * 2 ^ 75T "2 5" C23 

5. r?(fc) + N(p) -► <j(A/) + iV(p') 

1/(5) = K 5 + n 5 (C24) 



with 



K 5 = ig*NN^S N (p + k) ^75 (C25) 
H 5 = ig„NN^ ^ 5 S N ( P - fc') (C26) 



6. <r(fc) + N(p) <j(A/) + AT(p') 

t/(6) = K 6 + H 6 (C27) 

with 

K 6 = gl NN S N (p + k) (C28) 
= g 2 aNN S N (p - k') (C29) 

7. 7r (k,i)+N(p)^ p(k',j)+N(p') 

V(7) = Vl + V b 7 + K 7 + + VJ (C30) 



with 



K 7 = i—gpNNT^SNip + k) far* (C31) 

V b 7 = i f -^g pNN ^S N {p - k')F p , (C32) 

T - /7 /Wat i (g - fc) • <y #75 , . 

K = g^ijir ^ — (C33) 

m n q l — m% 

Vj = -^ 9p nn fylsejuT 1 (C34) 



m 



vj = aj^^g^ ^^hfy b , + ^_ ^ (C35) 

Ulu q — TYl^ 4771 ^ 
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where 

IV = y[£ + £^(f P > r £)] (C36) 
F(8) = K 8 + V* (C37) 

with 

K 8 = i f -^g P NNT P >S N {p + k) ^75 (C38) 

777^ 

? s = ihniL gpNN ^ l5 S N (p - fc')r> (C39) 

777^ 

9. <7(fe) + /Y(p) p(k',j) + N(p') 

V(9) = V! + V b 9 (C40) 

with 

K 9 = g P NNgaNNT p >S N (p + k) (C41) 

H 9 = ^iv^7VJv^(p-A; , )r p , (C42) 

10. p (k,i)+N(p)^ p'(k',j)+N(p') 

V(10) = V™ + V?° + V™ (C43) 

with 

Va° + H 10 = g 2 pNN [T p/ S N (p + k)T p + r,SW(p - k')Y pl ] (C44) 

where 

Tp = ^ P ~ £^(? P V- V (C45) 



K 10 = f-^yp f P >~ t, ^ir l (C46) 



54 



11. ir(k,i)+N(p) ^ir'(k',j)+A(p') 
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v. 



14 



v„ 



14 



■ 9 pAAfirN A 

-%— |e 



AlalPp 



\4p 



KpAA 

4mA 



{{ p V-V{ p ))TiSf{p-k')Tlk' 



(C58) 



(C59) 



15. TT(k,i)+A(p)^TT'(k',j)+A'{ P ') 



1/(15) = K i& + H i& + K i& + v? (C60) 



with 



^15 = [AiVA]2 e /* . ^'^(p + fc ) £A . fc(T i)t (C61) 



n 15 = [— ] a ^T^a-J^Cp + *)[<*],, ^75TA (C62) 



V d 15 = ^AAV^T^{(^ ^) + T^((^+ ^ *~ ^+ • tA 

qz _ m z ArriA 



The baryon propagators in Eqs.(C8)-(C64) are 



(C64) 



S N (p) = -r^— , (C65) 
p — ttin 

- + + ~ ^ • <-) 

Eq. ()C66|) is the simplest choice of many possible definitions of the A propagator. It is part 
of our phenomenology for this rather complex coupled-channel calculations. 

Although the expressions Eqs.(C8)-(C64) look like the usual Feynman amplitudes, the 
unitary transformation method defines definite procedures in evaluating the time component 
of each propagator. For each propagator, the vertex interactions associated with its ends 
define either a "virtual" process or a "real" process. The real process is the process that can 
occur in free space such as A — > ttN. The virtual processes, such as the ttN — > N, ttA — > A, 
and 7rA — > N transitions, are not allowed by the energy-momentum conservation. The 
consequences of the unitary transformation is the following. When both vertex interactions 
are 'virtual', the propagator is the average of the propagators calculated with two different 
momenta specified by the initial and final external momenta. For example, the propagator 
of V a of Eq. (C9), which corresponds to 7i(k)N(p) — > iV — > 7r(k')N(p'), should be evaluated 
by 

S N (p + k) - ^[S N (p + k) + S N (p' + k')} 
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I (E N (p) + Ej^jyf - 7 ■ (p + k) + m N 

{E N {p') + Ej^Yrf - 7 ■ (gj + g) + mjy 
(^(p') + K(fc')) 2 - (p' + fc') 2 - ^ 



(C67) 



One see clearly that the denominators of the above expression are independent of the collision 
energy E of scattering equation Eq. (|B3j) and finite in all real momentum region. This is the 
essence of the unitary transformation method in deriving the interactions from Lagrangians. 
When only one of the vertex interactions is 'real', the propagator is evaluated by using 
the momenta associated with the 'virtual' vertex. For example, the propagator of Vj 1 
Eq.(C51), which corresponds to ir(k)N(p) — > A — > n(k')A(p') is S A (p' + k'), not S A (p + k) 
or [S A \p' ' + k') + S A '(jp+k)]/2. The structure of its denominator is similar to that of Eq.(C67) 
and hence the resulting matrix elements are also independent of scattering energy E and 
finite in all momentum region. The terms which have one 'real' and one 'virtual' vertex 
interactions are V^, Vj, V^ 1 . Their corresponding intermediate momentum variables have 
been correctly specified. The average, such as that of Eq.(C67), must be used in all other 
terms of Eqs.(C8)-(C64). We note that there is no propagator in Eqs.(C8)-(C64) which is 
attached by two real processes such as ttN — > A — > nN. Such real processes are generated 
from IV of the Hamiltonian and included in the resonant term tf 4B M , B , of Eq.dHJ). 

APPENDIX D: MATRIX ELEMENTS OF VERTEX INTERACTIONS 

We need to have partial- wave matrix elements of vertex interactions T^ n N, hp-,^, and 
h a ^ wn to evaluate the self-energy S a with a = ttA, pN,aN of Eqs.()59 |l -(|HT |l . and the one- 
particle-exchange interactions Z^^ A , Z^ wA , and Z^j %A , illustrated in FigJH] 

In consistent with the normalizations defined by Eqs. (j48|) . we write the matrix element 
of the a(p a ) — > Pipp) +7(p 7 ) vertex interaction / a ,/? 7 as 



where j a is the spin of the particle a, l a is the relative orbital angular momentum of the 
pair (/?, 7), n a = [(/ a (j/3j 7 )s a )]j a ] denotes collectively all quantum numbers specifying the 
interacting (j3, 7) pair, The momenta are related by relativistic kinematics 



< Pa, jarrijjfa^lp 1 ' a q' ' a ; jpm ip j 1 m h >= 5{p a -p' a ) 

x < Jam- ja \l a s a mi a m Sa >< s a m Sa \j Ph m h m h > fn a W a )Yi amia {q' a ) , (Dl) 



allm's 




P'p + P' 7 , 

P'p + Pa(p' a ,p'/3,x)p' a , 



(D2) 

(D3) 
(D4) 



where x — p' a ■ p'p, and 
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The form factor g na (q' a ) f° r P ~ > 7X71 an d ° ~ > 7X71 are taken from Ref. [68j and A 
from SL model. They are related to phase shifts by 



-[sm5 a (E)]e i5 "(E) 



with 



E a (E) = Jq 2 



\fnM\' 



E-M°-V a (Ef) 

\fnM)? 



E — Ep(q) — E 7 (q) + ie 



Explicitly we have 



fA,TvN(q) 



NA 



E N (q) + m N A-k q 



'^{^f ,2 ^2EM)\ 2E N (q) V 3 ff[ <r» + A»* A J 
with f nNA = 2.049, A^va = 3.29 (fm) -1 , and M A = 1299.07 MeV, 

f ( \ - qr 

with g p ^ = 0.6684 m^ 2 , r = 0.428 fm, and M° = 811.7 MeV, 

1 

9a,7rn(q) = 9 a,; 



1 + (qr) 2 

with g a ^ = 0.7550my 2 , r = 0.522 fm, and M° = 896.8 MeV. 



ttN are 



(D5) 



(D6) 



(D7) 



(D8) 



(D9) 



APPENDIX E: MATRIX ELEMENTS OF Z^ >M , B ,(E) 

With the matrix elements of the vertex interactions defined by Eqs.(Dl)-(D4) of Appendix 
D, we can evaluate the partial- wave matrix elements of one-particle-exchange interaction 
%mb M'B' • We use the cyclic notation (a, (3, 7) to specify the particles involved in the vertex 
interaction a — > (3 + 7. By using the angular momentum quantum numbers defined in 
Appendix D, we then define the basis state of a MB system with a given total angular 
momentum (JM) in the center of mass frame as 

\N a ;p a JM >= \ {L a [(l a (j^)s a )] a j a }S a }JM;p a > , (El) 

where we have introduced a concise notation AL = [{L a [(l a (jpj 7 )s a )j a j a ]S a }]. 
Following the standard procedures of Ref.|69j], one then obtained 

< Np- P pJM\G^ N (E)\N a , Va JM > 

L a=0 6=0 
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where 8 a> p = 1 — 5 a ^, and 



(2l a + 1)1(2/^ + 1)! 1/2 
1 (2a)!(26)!(2/ Q -2a)!(2/ /3 -26)! t 



7l2 ) S a S3 f 



jot S a Sp j/3 



/AA' I L ^ L « J i [ S Q Z a S/3 



X 



^? ^ { W a i fl -6 A 



l a — a b A' 

, A' L L p \ ( A L L Q \ / Z a - a b A' \ I a lp - b A . , 
l 000 000 00 o o ' *' 



with a = \/2a + 1 and 

R = -J + L + L a + Lp + S a + Sp + j a + jp + sp + l p - j a . (E4) 

In the above equations, the usual 3-j, 6-j, 9-j, and 12-j symboles have been used to define 
the angular momentum coupling. The details can be found in Ref.j6^. 

The three-body cut effects are in F£ na (pp,p a ] E) of Eq.(E2). They are calculated from 
the vertex functions f na (q a ) by 

LAa . , 1 r« Qp l "pp(x)fl( q p) q -^pUx)f na (qa)PL(x) 
t n l) n a \Vp,VwE) - - / dX— T—i / — » \ r-i / -> , n , • , (E5J 

2J-i E- E a ,{-p a ) - Ep,{-pp) - E y (p a +pp) + te 

where x — p a ■ pp and the vertex function f Ha defines a — > [3' + 7 and f n „ defines 7 + a' — > (3. 
Namely, a' and (3' are the spectators of the decay of particle a and (3 respectively. We 
obviously have a = A, a' = ir, 7 = N, (3 = A and /3 — w for Z^ A (E), and a = p, a, 

a' = N, 7 = 71, (3 = A and (3' = vr for ^ wA (£). 

In the actual calculations, the integration path — 1 < x < 1 of eq.(E5) is deformed into the 
complex x— plane in order to avoid the singularity x$ (—1 < Xq < 1) where the denominator 
vanishes. We have used a simple parabolic form, i.e., x = t + i(t 2 — 1). 



APPENDIX F: MATRIX ELEMENTS OF -yN -» MB TRANSITIONS 

To include the final meson-baryon interactions in the photo-production, it is only neces- 
sary to perform the partial-wave decomposition of the final MB state. We thus introduce 
the following helicity-LSJ mixed-representation 



(k', q)= Ely < ^V' M (-A' B ) I S'S' Z > 



V L'S'M'B',X 1 X N 



x 'm x 'b 



x < L'S'0S' z I JS' Z > 

x < J, k'X' M (—X' B ) I VM'B'nN I J,q\{-X N ) >] , (Fl) 



59 



where < J,k'X' M (—X' B ) \ vm'B'^n \ J, l^-yi— \n) > can be evaluated using the same expres- 
sion of Eq.(C2) using the heliclity matrix elements of vm'B'^n- To evaluate these quantities 
with the normalization defined by Eq. (j48J) . we define for a photon four momentum q = (u},q) 

< (k'j),p'\v MB „ N \q,p> = -}= < (k'j),p'\J2 ^ M (™)e M |<?,P > 

V Z( l0 n 



{2-Kf^^E' B {k')\2E' M {k') ^ ' ^'\jE N {q)^2q- Q 

(F2) 

where is the photon polarization vector, and n denotes a given considered process 

I(n)=e-](n). (F3) 

Here j(n) can be constructed by using the Feynman rules. The resulting expressions for 
each of 7iV — > ttN, rjN, aN, pN, 7rA are listed below : 

1. 7 (q) + N(p)^Tr(k',j)+N(p') 



with 



-f(l) = 'l + 'I + 'l + 'I +'} + 'l + 'l (F4) 



where T N = e N ^ - -^-[^ 4~ 4 ^y] (F6) 

il = ii!™r N —± JfV (FT) 

m n p— f — m N 



Il = i^^Tt m '^ST(p-k')k'T j 



Il = ^e tj3 r* (F9) 

771,1- 

I e = j- -e ij3 T (k + k ) • e 7 (F10) 

where k = p — p 

xe^'kVe^j-L- (Fll) 
k A — nip 



(F8) 
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! = _ 9 u NN9u»n [ s + J^_ { S fogy 

y m w Ann jv 

xe a p vS k ri q a e l ?5 j3 ~ 



k 2 — m 2 



where T e ™ A = T e ™' A e». 



2. j(q) + N(p) 17 ( A;') + TV(p') 



J(2) = / a 2 + / 6 2 + J 



with 



j2 _ JyNN p 1 

ft! — 

il = i^Lr N —l ft' l5 

/2 = _ 9pNN9 f rn f_ [ u+ _Kp_, „ u_ u v^ 

m p 2 Am N 



3. 7 (g) + N[p) a(fe') + iV(p') 



with 



A3) = 4 3 + / b 3 



p'+ ft' — m N 

lb — —QaNN^N- 



ft' — itln 

4. 1 ( q ) + N( P )^p(k',j,X) + N(p') 
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-g P N N r p/ — p_ m /N (F21) 

-^N N T N ^_^_ m Y pl (F22) 

^(K f P >- VfaH^SZW + k')YT A (F23) 

Ulp 

_ fpNA [rr A ]t ^ (p _ k r )T J {K f _ re * ;j7s (p24) 



< 1^ C <W£ ^- ^ 7 £] (F25) 

x + k>) ■e.-ik- e* p ,)e» - (e 7 • f)k'*]J^ ( F 26) 

= -^^V ^^fc'^Vc?^-^ (F27) 

(F28) 



7r(fc',j)+A(p') 



7(5) = If + J, 5 + J c 5 + 7| + J e 5 + I) + I 5 g (F29) 



il^ e * A ■ k'T^S N (p' + k')T N (F30) 

i UNA r emA e ^ SN{p _ k >) ^ ^ (F31) 

-i—^ fatfiSZtf + k')Tl m > A (F32) 

.^rA^ ( 1 + T 3 )h< ^ + (£y)fj7l5 A (p _ (F g 3) 

I I Iff £ 

^^sT\ ■ el (F34) 

/vV %3r[lf + ^] (F35) 
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Il = -^^^f ■ *f ~ ~ k " A\l^<?e^ (F36) 

where the pion pole term Ig consists of energy independent interaction V g 5 and energy 
dependent interaction Z\ given as 



1 e* A ■ ki(ki + k') ■ e 7 



9 ~ 2E n {k - k') E N {q) - E A {k') - E n {k - k') A 7 1 ' 



Z 5 = 1 4 • fe 2(fc 2 + k') ■ 6 7 

9 2E n (k-k')E-E N (q)-E n (k')-E w (k-k')+ie ' 



with fci = {E n (k — A;'), fc' — fc) and &2 = (—E n (k — k'), k' — k). The on-shell matrix element 
of V^ 5 + Zg is given as 

V 5 + Z 5 = e* A ■ k(k + k') ■ e 7 1 (F39) 

«r — ml 



APPENDIX G: MULTIPOLE AMPLITUDES OF 7 iV -» vriV 

For 7iV — > Mi? matrix elements, we use the helicity-LS J mixed-representation defined by 
Eq.fjBU)- It can be calculated by using Eq.(Fl). For pseudo-scalar meson 7r and rj production, 
it is often to write the amplitudes in terms of multipole amplitudes. Here we want to relate 
our matrix element Eq.(Fl) and hence also the amplitude Eq. (j64J) to this commonly used 
multipole amplitude. 

With the definition Eq. (|TT|) for the scattering amplitude T, we first define the amplitude 
F by the on-shell T-matrix element of 'yN — > MB as 

< M5|F| 7 iV > = -^. y /E N (k)E M (k)\q \E N (q)^L= < MB\T\ 7 N > (Gl) 

w V2M 

with 

T = J% = J°e° - J ■ e, (G2) 

where J M e M = J2 n j^(. n ) e n is identical to that in Eq.(F2). The most general Chew- 
Goldberger-Low-Nambu (CGLN) amplitudes [TtJ F can be written as (isospin index is sup- 
pressed) 

F = —icr ■ e±F\ — cr ■ ker ■ q x e\_Fi — icr ■ qk ■ €±F% — icr ■ kk ■ e±F^ 

—icr ■ qq ■ eF 5 — icr ■ kq ■ eF 6 + icr ■ ke F 7 + icr ■ qe F 8 . (G3) 

Each coefficient in the above equation can be written in terms of multipole amplitudes E\±, 
Mi±, Li± } Si± 

F x = J2[P! + xEi + + PUEi- + lPl +x M l+ + (I + l)PUM^] , (G4) 
i 

F2 = £[(/ + X ) p / M '+ + lP ! M i-] > (G5) 
i 
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C 1 

-T3 


\^rp" j? j_ p" p p" A/f 

1 


_L P" A/f 1 




771 


1 


P" A/T 1 
M l-\ > 
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(G8J 


F 6 






(G9) 






(G10) 








(Gil) 



where Pl{x) is Legendre function and x = k ■ q. It is of course well known that only four of 
the above amplitudes are independent for photo-production and six for electro-production. 
Choosing the photon direction as q — z, Eqs.(F2) and Eqs.(Gl)-(G2) clearly lead to 



rrr-r. ttt—tt-t 1 



<J|F|A 7 s> = - — ^E N (k)iu M (k)\q \E N (q)^—vi^ x Jk,q)^2q . (G12) 



Here s is the z-component of the initial nucleon spin and we have dropped the notation 
MB = irN and (LS) — (I — J ± 1/2, 1/2) and isospin T in defining the matrix element of 
F. 

With the form Eq.(G3), it is easy to calculate the matrix element < J|F|A 7 s > in our 
helicity-LSJ mixed-representation. After some derivations, we obtain the following relations 

E i+ =5sfelj[ <J = l + l/2\F\X=l, 8 = -l/2> 



< J = l + l/2\F\\= l,s = 1/2 >] (G13) 



V 1 + 2 

Ei- =jtul -<J = l-l/2\F\X=l,s = -l/2> 



? ' 1 <J = l- 1/2\F\X = 1, s = 1/2 >] (G14) 



V 1-1 

M '+ =4^7TT)[ <J = l + l/2\F\\=l,s = -l/2> 



+ \J l ~j^< J = l + l/2\F\\ = l,s = l/2 >] (G15) 
Mi- =^u\ <J = l-l/2\F\\ = l,s = -l/2> 

< J = l-1/2\F\\ = l,s = l/2 >] (G16) 

^+ = -i^TT) < J = 1 + V2|F|A = 0, 8 = 1/2 > (G17) 

L,_ = < J = Z-l/2|F|A = 0,s = l/2> (G18) 

^+ =4^TT) <^ = / + l/2|F|A = t,s = l/2> (G19) 

Si- = <J = l-l/2\F\X = t,s = l/2> (G20) 
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Here we used polarization vector e M (A = t) — (1, 0). Substituting Eq.(G12) into Eqs.(G13)- 
(20), we can relate the usual multipole amplitudes of 7iV — > ttN to the matrix element 
Eq.(Fl) in the helicity-LSJ mixed-representation. 



APPENDIX H: jN vrvriV AMPLITUDES 

We consider the non-resonant ~f(q)N(p) — > Tr(k l )7r(y)N(p') illustrated in FigEJ With 
the normalization defined by Eq.(|73p. we define the matrix element of this amplitude with 
a photon momentum g M = (u,q) as 

< k j (j),k\i),p | twv, 7 7v | p >= < I J I P > -e 7 (?) > ( H1 ) 

V2u; 



where z and j denote the isospin components of the produced nn, e 7 is the photon polariza- 
tion vector, and 



< W(j),k l (i),p' | J M | p > 



(2vr)9/ 2 V ^(P') ^JJ^EjkT)\ E N (p) 
x^(k')j^(k) (H2) 

with 

= + f{2) + ^(3) + j"(4) + ^(5) + j"(6) (H3) 

Each term of Eq.(H3) are from mechanisms illustrated in Fig0] Within our formulation, 
the non-resonant mechanisms are only from diagrams with intermediate nucleon states. The 
exchange mesons can be tt, p and u. We then have the following expressions 

= i{ Unn^ + fc <)y* 7se r * + ^ l5 e kj3 r k S N (p - fc*) tfV*] , (H4) 

+ Pl^SNip' + k l )r N S N {p - k j ) ft l5 r j 

+J^S N (p — k i — tf) f( l l5 T l S N (p - k j ) , (H5) 

+{*- /" F)l 5 e kj3 r k S N (p - k l ) ^ ^I^r-Zl ] ' (H6) 



m = ^i UNN9p T 9p ^ ]{^s N ( P > + k^tf + -^( 7 3 ?- ^)] 

m% 2 4m jy 

r J ' r x K n , x r. 7, X, „ , , „-s .„• 1 



+ -[7 d + T^(7 d #7 d )SW(p - k l ) r7 5 r 1 ]]^ , (H7) 

l 2 



2 4mAr P - m* 



m_ 4mjv 



4m 
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Xe a p vS Wq a e% 3 T 1 , (H8) 
k z — 

m = -g P N N g P J lh * + t36 ; 3 - 2T3kj ir + ^mv i- h% — L — 2-( H9 ) 

2 4 771 7v (p — p) z — m z l 

In the above equations, we have defined 

r N = l A^r+i^^^ qv , (hio) 

and 

k—p — p' — k 1 . (Hll) 



APPENDIX I: RESONANT AMPLITUDES 

In this appendix, we give formula for using the N* parameters listed H2] by Particle Data 
Group to calculate the resonant amplitudes defined by Eq.([BT)|). 

In the rest frame of iV*, the amplitudes of strong decays of a N* with a mass M JT and 
spin-iospin (JT) can be written as the following partial- wave form 

< 4>j mj ,Tm T | T N *^ MB | k,m jM m tM ,m jB m tB > 

= EE[ < Tm T \t M t B m tM m tB >< Jmj\LSm L m s >< SM s \j M j B m jM m jB > 

LS all m z 



where k B is defined by M JT = E B (kn) + E^(kji) an d the form factor is chosen such that 
fis(kn, fcjj) = 1. With the normalizations 

< fc|A? >= 5(A? - y?) , (12) 
the partial decay widths can be written as 

dT MB (N* JT ) = (2ti)5(M jt - E B (k) - E M {k))^-dk 

E l< 4>Jmj,TmT I T N *^ MB | k,m jM m tM ,m jB m tB > | 2 ] . (13) 

mj m jM ,m jB 

From Eqs. flllt and (|j*3)). we then have 

TM B (iV} T ) = E|Gi|| 2 . (14) 

Eq. ffl4|) allows us to determine the coupling constant up to its phase in terms of the 
empirical partial decay widths as listed by Particle Data Group [62|. Here we use the phase 
from the 3 P model of Capstick and Roberts [Til]. 
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For N* — > *yN amplitudes, we use the commonly used helicity representation to define 
dT yN (N} T ) = (2v)5(M JT - E N (q) - q)^-^dq 

EE l< 

Jmj ,Tmj< I Fn*->mb I q, A 7 , Aat, m 4jv > | 2 ] (15) 

"V A 7 Ajv 

With the normalizations defined by Eq. (jl2|) . we then define 

= ^ w tf^_ Aw) ^ (16) 

where g\ r {q,qR) is a form factor with g# defined by M JT = qn + E^{q^) and normalized 
as 9\ T {qR,qR) — 1- Substituting Eq. (Jl6j) into Eq. (Jl5j) and noting that \A J _J X \ = A J X T , we then 
obtain the standard form 

We only include 3 and 4 stars iV* in our calculations of Eq (j66p . We use their mean values 
of G J L T S and A x , as listed in Tables EKD 
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TABLE IV: The helicity amplitude A x is given in unit of 1(T 3 GeV^ 1 / 2 . G LS is in unit of MeV 1 / 2 . 
The resonance mass and the total decay width T tot are in unit of MeV. 



N$j{M R ) 




channels 


L,5 


Gls 


A/2 


^3/2 


Si i f 1535) 


150 








0.090 


0.0 






ttN 


1/2 


6.26 










nN 


1/2 


7.55 










7rA 


2 3/2 


1.06 










pN 


1/2 


1.49 










aN 


1,1/2 


1.50 


- 


- 


Sn (1650) 


150 








0.063 


0.0 






ttN 


1/2 


12.23 










nN 


1/2 


3.48 










7rA 


2,3/2 


2.01 


- 


- 






pN 


0, 1/2 


1.42 












2,1/2 


5.124 










aN 


1,1/2 


1.42 






Pi (1440) 

■*■ 11 V v / 


350 


^N 






-0.065 


0.0 






ttN 


1 1/2 


18.78 










7rA 


1,3/2 


8.85 


- 


- 






aN 


0, 1/2 


7.66 






Pi i f 1710) 


100 


^N 






0.009 


0.0 






7riV 


1, 1/2 


6.22 










7]N 


1,1/2 


2.93 










7rA 


1,3/2 


7.47 












1,1/2 


4.93 










aN 


0,1/2 


1.19 






Pi 3 (1720) 


150 


7 iV 






0.018 


-0.019 






vriV 


1,1/2 


2.45 












1,1/2 


2.20 







pN 1,1/2 10.49 
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TABLE V: The helicity amplitude A x is given in unit of 1CT 3 GeV^ 1 / 2 . G LS is in unit of MeV 1 / 2 . 
The resonance mass and the total decay width T tot are in unit of MeV. 



N}j(M R ) 


ptot 


channels 


L,S 


Gls 


A/2 


As/2 




Di 3 (1520) 


120 








-0.024 


0.166 






ttN 


2, 1/2 


8.84 










ttA 


0,3/2 


4.31 












2,3/2 


3.69 










P N 


0,3/2 


3.34 










aN 


1,1/2 


1.11 






013(1700) 


100 


7 iV 






-0.018 


-0.002 






7TN 


2,1/2 


2.65 










ttA 


0,3/2 


4.38 












2,3/2 


11.758 










pN 


0,3/2 


3.5 








015(1675) 


150 


7 iV 






0.019 


0.015 






7TN 


2,1/2 


6.77 










ttA 


2,3/2 


9.085 










pN 


2,3/2 


1.46 








F 15 (1700) 


130 


7 iV 






-0.015 


0.133 






TTN 


3,1/2 


9.39 










7rA 


1,3/2 


4.23 












3,3/2 


1.13 










pN 


1,3/2 


2.52 












3,3/2 


1.95 










aN 


2,1/2 


3.39 






Gi 7 (2190) 


450 


7 iV 






-0.055 


0.081 






7TN 


4,1/2 


9.52 










pN 


2,3/2 


11.46 
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TABLE VI: The helicity amplitude A x is given in unit of 1(T 3 GeV^ 1 / 2 . G LS is in unit of MeV 1 / 2 . 
The resonance mass and the total decay width T tot are in unit of MeV. 



N}j(M R ) 




channels 


L,S 


Gls 


A/2 


A3/2 




5 3 i(1620) 


150 


7 iV 






0.027 








irN 


0,1/2 


8.02 










vrA 


2,3/2 


7.47 










pN 


0,1/2 


4.57 












2,3/2 


1.69 








^31 (1910) 


150 


7 iV 






0.003 








ttN 


1,1/2 


14.38 










1,1/2 


11.5 










P 33 (1600) 


350 


7 iV 






-0.023 


-0.009 






ttN 


1,1/2 


11.75 










vrA 


1,3/2 


17.06 








P 33 (1920) 


200 


7 iV 






0.04 


0.023 






vriV 


1,1/2 


2.48 










vrA 


1,3/2 


7.10 








P 33 (1700) 


300 


7 iV 






0.104 


0.085 






irN 


2,1/2 


2.44 










vrA 


0,3/2 


10.35 












2,3/2 


2.18 










piV 


0,3/2 


1.09 








P 35 (1905) 


350 


7 iV 






0.026 


-0.045 






vriV 


3,1/2 


6.11 










vrA 


1,3/2 


9.78 












3,3/2 


13.53 










pN 


1,3/2 


9.99 








P 35 (1950) 


300 


7 iV 






-0.076 


-0.097 






TTiV 


3,1/2 


10.38 










vrA 


3,3/2 


9.39 
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